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Preface 


Hydraulics is a branch of scientific and engineering discipline that deals with the mechanical properties 
of fluids, mainly water. It is widely applied in many civil and environmental engineering systems 
(water resources management, flood defence, harbour and port, bridge, building, environment 
protection, hydropower, irrigation, ecosystem, etc). This is an introductory book on hydraulics and 
written for undergraduate students in civil and environmental engineering, environmental science and 
geography. The aim of this book is to provide a concise and comprehensive coverage of hydraulics that 
is easy to access through the Internet. 


The book content covers the fundamental theories (continuity, energy and momentum equations), 
hydrostatics, pipe flow, physical modelling (dimensional analysis and similarity), open channel flow, 
uniform flow, channel design, critical flow, rapidly varied flow, hydraulic jump, hydraulic structures, 
gradually varied flow, computation of flow profile, unsteady flow and hydraulic machinery (pump and 
turbine). The text has been written in a concise format that is integrated with the relevant graphics. 
There are many examples to further explain the theories introduced. The questions at the end of each 
chapter are accompanied by the corresponding answers and full solutions. A list of recommended 
reading resources is provided in the appendix for readers to further explore the interested hydraulics 
topics. 


Due to its online format, it is expected that the book will be updated regularly. If you find any errors 
and inaccuracies in the book, you are encouraged to email me with feedback and suggestions for 
further improvements. 


Dawei Han 
Reader in Civil and Environmental Engineering, 


Water and Environmental Management Research Centre 
Department of Civil Engineering 

University of Bristol, BS8 1TR, UK 

E-mail: d.han@bristol.ac.uk 
http://www.bris.ac.uk/civilengineering/person/d.han.html 


August 2008 
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1. Fundamentals 


Hydraulics is a branch of scientific and engineering discipline that deals with the mechanical properties 
of fluids, mainly water. It is widely applied in water resources, harbour and port, bridge, building, 
environment, hydropower, pumps, turbines, etc. 


1.1 Properties of Fluids 


1) Density: Mass per unit volume p (kg/m”) 
For water, p =1000kg/m* at 4°C, 
p =998 kg/m? at 20°C 
For air, p=1.2kg/m? at 20°C at standard pressure 


2) Specific gravity: Ratio of the substance’s density and water’s density at 4°C 


3) Pressure: Normal fluid force divided by area over which it acts (N/m2). (Note: pressure is scalar 
while force is a vector. A force is generated by the action of pressure on a surface and its direction is 
given by the surface orientation) 


4) Viscosity and shear stress 


Take an element from the fluid 


Moving with U cae 
Stationary 
du 
T= UU 
dy 


and the total force 


UA 
[ne ea 
a 


where p absolute viscosity (N s/m’),7t shear stress (N/m”). Fluids that follow the aforementioned 
formulas are called Newtonian fluids. 
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Examples 

Air ut =1.79 x 10° Ns/m?* 
Water w= 1.137x10° Ns/m? 
Engine oil SAE30W u=0.44 Ns/m* 


High viscosity: sticky fluid; low viscosity: slippy fluid 


u is not a constant and changes with temperature. 


B= 


Another form of viscosity is kinematic viscosity v= (m’/s) 


p 
Assumptions for the equation: 


1) Fluids are Newtonian fluids (Non-Newtonian fluids are studied by Rheology, the science of 
deformation and flow); 


2) The continuum approximation: the properties of the fluid can be represented by continuous fields 
representing averages over many molecules (The exception is when we are dealing with gases at low 
pressures). 


1.2 Flow Description 


There are two approaches 

1) Lagrangian approach: follow individual fluid element as it moves about; 

2) Eulerian approach: focus on a fixed location and consider how the fluid properties change at that 
location as time goes on. 


Definitions relating to Fluids in Motion 
Ideal flow: frictionless and incompressible (i.e. nonviscous). 
Steady flow: The flow is steady if the properties at each point in the fluid do not change with time. 


One, Two and 3D flows: 

One dimensional flow requires only one coordinate to describe the change in flow properties. 
Two dimensional flow requires two coordinates to describe the change in flow properties. 
Three dimensional flow requires all three coordinates to describe the change in flow properties. 
In general, most flow fields are three dimensional. However, many practical problems can be 
simplified into one or two dimensions for computational convenience. 


Download free ebooks at bookboon.com 


10 


Please click the advert 


Concise Hydraulics Fundamentals 


Control Volume 


A suitable mass of fluid can be identified by using control volume. A control volume is a purely 
imaginary region within a body of flowing fluid. The region is usually at a fixed location and of fixed 
size. Inside the region, all the dynamic forces cancel each other. Attention may therefore be focused 


on the forces acting externally on the control volume. 
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1.3 Fundamental Laws of Physics 


The fundamental equations that govern the motions of fluids such as water are derived from the basic 
laws of physics, i.e., the conservation laws of mass, momentum and energy. The conservation of 
momentum comes from Newton’s second law of motion stated in 1687 in Principia Mathematica. The 
law of conservation of mass was formulated in the late eighteen century and the law of conservation of 
energy in the mid-nineteenth century. In modern physics, mass and energy can be converted from one 
to the other as suggested by Albert Einstein in 1905. This combines two individual conservation laws 
into a single law of conservation of mass/energy. However, since conversion between mass and energy 
are not of relevance to fluids studied by hydraulics, the two individual laws of conservation are used in 
hydraulics. The conservation of mass is also called the equation of continuity. 


Conservation of mass: mass can be neither created nor destroyed. 

Conservation of energy: energy can be neither created nor destroyed. 

Conservation of momentum: a body in motion cannot gain or lose momentum unless some external 
force is applied. 


The application of the three fundamental laws will be further explained in the following chapters. 


Questions 1 
Fundamentals 
1. A plate separated by 0.5 mm from a fixed plate moves at 0.5 m/s under a shear stress of 4.0 N/m’. 
Determine the viscosity of the fluid between the plates. 


(Answer: 0.004 Ns/m’) 


2. A Newtonian fluid fills the gap between a shaft and a concentric sleeve. When a force of 788N is 
applied to the sleeve parallel to the shaft, the sleeve attains a speed of 2m/s. If a 1400 N force is 
applied, what speed will the sleeve attain? The temperature of the sleeve remains constant. 


(Answer: 3.55 m/s) 


3. Water is moving through a pipe. The velocity profile at some section is shown below and is given 
2 
2 
mathematically as U = Awl 4. —T |, where u = velocity of water at any position r, B = a constant, 
L 


4 


= viscosity of water, d = pipe diameter, and r = radial distance from centreline. What is the shear stress 
at the wall of the pipe due to water? What is the shear stress at a position r = d/4? If the given profile 
persists a distance L along the pipe, what drag is induced on the pipe by the water in the direction of 
flow over this distance? 
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(Answer: d/4, d/8, Bd? 1L/4) 


Solutions 1 

Fundamentals 
1. F a 
. From t= u— 
Ca 


7 0.5 
~ 0.5 x 0.001’ 
so u = 0.004 Ns/m2 


4.0 


5 F U 
.T= =U 
Ad 
FA 
R — = — =constant 
ealrange as Ud 
Fi_#, 
Therefore, = 
Uo 
788 _ 1400 
a” ae 
So U, =3.55m/s 
2 
3. From u = at - 2) 
4ui4 J 
so La a ee ee 
dr 4u 2 
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Then 
At the wall, r= d/2 , 


hence 7,7 a 


? 


Tad 14 = 


Drag = T,, (area) -= (dL) = fd’ a 14 


<< 
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2. Hydrostatics 


Hydrostatics deals with the pressures and forces resulting from the weight of fluids at rest. 


2.1 Pressure 


Pressure is the force per unit area acting on a surface and it acts equally in all directions. A common 
pressure unit is N/m’, also called Pascal (Pa). Depending on the benchmark used (with/without 
atmospheric pressure), pressure can be described as absolute pressure or relative pressure. 


1) Atmospheric pressure ( p, ): is the pressure at any given point in the Earth's atmosphere caused 


by the weight of air above the measurement point. The standard atmosphere (symbol: atm) is a 
unit of pressure equal to 101.325 kPa (or 760 mmHg, 1013.25 millibars). 


2) Absolute pressure ( p,,,, ): is the pressure with its zero point set at the vacuum pressure. 
3) Relative pressure ( p, ): is the pressure with its zero set at the atmospheric pressure. It is more 


widely used in engineering than absolute pressure. 


The relationship between them is 

P, = Pars ~ Pa 

The change of pressure within the fluid can be expressed as 
dp = —pgdz 


For a fluid with constant density, the differential formula can 
be integrated as 


z+ =€ 


Pg 


where C is a constant that can be set from the boundary 
condition. If the top boundary pressure is py and depth of the 
fluid is h, the bottom pressure can be derived as 


p=pgh+ p, 
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For any two points in the same fluid, it can be derived 


z,+24 a7, 4 De 


pg 


2.2 Manometer 


A simple manometer is a tube with its one end attached 

to the fluid and the other one open to the atmosphere = 
(also called Piezometer). The pressure at Point A can be 
derived from the height h, in the tube. 


Pa = pgh, 
A more complicated manometer uses U- 


tube with a different fluid to the measured 
one. The pressure at Point A will be 


Pa = PmGhy — PGh, 
U-tube manometer Differential manometer 

A differential manometer is used to measure the pressure difference between two points. A U-tube with 

mercury (or other fluids) is attached to two points whose pressures are to be measured. The pressure 

difference can be derived by measuring the elevation differences between Point A and Point B (i.e., Az), 

and between the mercury levels (Ah). 


Therefore 
Pa Pg = Pn — p)Ah + pAz] 


2.3 Pressure Force on Plane Surface 


For a plane surface with area A, the total pressure force can be derived by the following integration 
formula: 


F =| pdA=| pghdA 
A A 


where h is the depth of fluid from its surface. 


If the centroid of the area is known, the 
pressure force can be derived as 
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F = pghA 

where fh is the depth of A’s centroid. 

The centre of pressure is the point through which the resultant pressure acts. 
[ypdA [y’aa 
A A 


- [paa [yaa 
A A 


y 


The numerator is the moment of inertia of the surface about the axis through O, and it equals 


tol, =I, +y-A. Therefore 
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_yyte 
sey YA 
Shape Area A | Centroid y | Moment of Inertia I, 
Circle (r radius) ar? r ar’ /4 
Rectangle (b width, h height ) bh h/2 bh? /12 
Triangle (b bottom width, h 
height) bh/2 2h/3 bh? /36 


2.4 Pressure Force on Curved Surface 


For curved surfaces, the pressure force is divided into horizontal and vertical components. The vertical 
force F', is the total weight of the fluid above the curved surface and its centre of pressure acts through 
its centre of gravity. The horizontal force F, equals to the pressure force on a vertical plane surface 
projected by the curved surface. The resultant force is a triangular combination of the horizontal and 


vertical parts. 


F,=pgV_ where V is the volume of the fluid above the curved surface. 


F=/F?+F? 


2.5 Flotation 
1) Buoyancy 


Buoyancy is the upward force on an object from the surrounding fluid (liquid or gas). This force is due 
to the pressure difference of the fluid at different parts of the object. The buoyancy force is derived by 
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F,=pgV 
where V is the displaced fluid volume by the object, p is the fluid density. 


An object submerged in a fluid is subject to two forces: gravity and buoyancy. When gravity force is 
greater than buoyancy force, the object sinks to the bottom. If the situation is reversed, the object will 
float. If two forces equal to each other, the object could be anywhere in the fluid. 


2) Stable flotation and metacentre 


A floating object is stable if it tends to return to its original position after an angular displacement. This 
can be illustrated by the following example. When a vessel is tipped, the centre of buoyancy moves 
from C to C’. This is because the volume of displaced water at the left of G has been decreased while 
the volume of displaced water to the right is increased. The centre of buoyancy, being at the centre of 
gravity of the displaced water, moves to point C’, and a vertical line through this point passes G and 
intersects the original vertical at M. The distance Gm is known as the metacentric height. This 
illustrates the fundamental law of stability. When M is above G, the metacentric height is positive and 
the floating body is stable, otherwise it is unstable. 


we 
I 


| 
| 
| 
| 
i 


iu 


I 


GM sina 


There are two ways to find the metacentre: experiment and analytical method. 


a) Metacentre by experiment (if a ship is already built, the experiment method is easy to apply) 
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Shift a known weight o from the centre of the ship by a distance / to create a turning moment P=oal 
and the ship (with a total weight of W) is tilted by an angle a. The metacentre can be derived by the 
balance of moments at the point G. 


GM = P = ol ey OO l 
Wsinag Wsina Wa 


b) Metacentre by analytical method 


Image the displacement centre (centroid of the buried body) is 
moved by x due to a tuning moment. This centroid displacement is 


contributed from only the top two triangles (worked out one triangle 
and the other one is doing the same thing, with either added buoyancy and reduced buoyancy). If the 
ship’s displacement volume is V, length is L and width is D, we can derive around C point the 
following, 
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D/2 
Vx = 2 } L(xtan @)xdx =(tana)I,, (1, is the 2" moment of the area) 
0 


ae (tana)I, 
Vv 


Therefore 


From the centre of gravity and buoyancy, 


GM =CM —-CG 


Questions 2 

Hydrostatics 
1. A rectangular plate gate is placed in a water channel (density of water: 1000kg/m’). Its width is 0.8m 
and the water depth is 2m. Estimate the pressure force and its centre of pressure 


(Answer: 15.7 kN, 1.33 m) 
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2. Estimate the resultant force of water on a quadrant gate. Its radius is 1m and width 2m. The centre 
of gravity for the quarter circular sector is 4R/3z (from the circle centre to the right). 


(Answer: 18.3 kN, 57.3°) 


3. A cube of timber 1m on each side floats in water. The density of the timber is 600 kg/m’. Estimate 
the submerged depth of the cube. 


(Answer: 0.6 m) 


4. A crane barge which is rectangular in plan, has a large water plane area compared with its depth for 
stability. The barge weighs 500 kN and the crane 50 kN. The resulting centre of gravity of the 
combination is at the level of the deck, while the centre of gravity of the crane itself is 3m above the 
deck. 


Tests are being undertaken to ensure the stability of the crane barge. The crane is moved horizontally 
sideways by 0.8 m and the barge rolls through an angle of 5°. What is the metacentric height of the 
system? When the crane is back in its central position, we need to know how high the jib can be raised 
before the barge becomes unstable. 


(Answer: 0.834 m, 9.2m) 


Solutions 2 
Hydrostatics 


1. The centroid of the gate is 


The area is 
A=bh=0.8x2=1.6 m* 


Therefore 
F = pghA=1000x9.81x1x1.6 =15.7 kN 


The centre of pressure 
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_ bh? _ 0.8x 2 


o=—= = 0.533 m* 
12. 12 


The gate is vertical, so y and h are the same. 


h= ee ce ee 
hA 1x1.6 
2. a) Vertical force 
mx 


F, = pgV =1000x9.81x 


Its centre of pressure is at 
4R/3 =0.424 m to the right of the circle 


b) Horizontal force 


x2=15.4 kN 


centre. 


The centroid for a vertical plane surface of 1m tall and 2 m wide is 0.5m 


F. = pghA=1000x9.81x 0.5x 2 =9.81 kN 
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Its centre of pressure is at 1/3 of the height 
0.333m to the bottom. 


c) Resultant force 


F =,/F?+F? = 15.4? +9.81 =18.3 kN 
Direction 
6 = tan '(F, / F,) = tan™'(15.4/9.81) =1.00(radian) = 57.3° 


It is 57.3° below the horizontal line and passes through the joint point between the vertical force line 
and horizontal force line.. 


0.424m 


3. Weight of the cube 
W = Py o0aI Vvood = 900 x 9.81 = 5886.N 


The volume of water to be displaced is 
WwW 5886 
Veta ~ = =0 6m° 


~ Poaeg (10009 


So the submerged depth will be 


= Veer 98 _ 9 6m 


Area 1x1 


4. This is to find the metacentre by experiment, so 
GM = S _ 50000 x 0.8 _ 0.834m 
Wsina 550000sin(5°) 


The metacentre is 0.834m above the deck. 


The centre of gravity of the barge is 
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500h = 3x 50 
h =0.3m (below the deck surface) 


If the centre of jib gravity is up by H from the deck and the system 


centre of gravity moves by 0.834m above the deck 
50 x (H — 0.834) = 500 x (0.3 + 0.834) 
H=12.2m 


So it can move by 
12.2 -3 =9.2m 


Plug into The Power of KnowledgeEngineering. 
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3. Energy Equation 


The energy equation is based on the conservation of energy, i.e., energy can be neither created nor 
destroyed. 


3.1 Basic Formula 
1) From solids 


Total energy = kinetic energy + potential energy (unit: joule) 


2 
Total energy = = 


+mqgz 


2) For fluids 
In comparison with the solid mechanics, there is an extra energy term in fluids (pressure). 


Total energy = kinetic + potential + pressure (unit: metre) 


. ue P sick 
Le. — te (per unit weight) 
28 8 


Pressure energy is similar to potential energy and they are closely linked (z+ p/ pg =C). One thing to 


be noticed is that in physics, energy unit is joule, but in hydraulics, engineers use meter and call it 
energy head. Unlike solid objects, fluid can move around and change its shape, so we use metre per 
unit fluid weight to describe the energy in fluid. One meter of energy head is equivalent of one joule 
energy of 1 Newton weight of fluid. 


3) Energy equation and continuity equation 


With the conservation of energy, the energy equation for one dimensional flow can be derived as 
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This is called the Bernoulli equation and it has been widely used in practical problems. 


In addition to the energy equation, the conservation of mass is usually used jointly to solve fluid 
problems. 


Rate of mass flow across 1 = Rate of mass flow across 2 
PAM, = pP,AM, 

If p is constant 
A, = Aw, 

It is usually called the continuity equation. 

3.2 Applications 

1) Example 

A large tank with a well-rounded, small opening 2 

as an outlet. What is the velocity of a jet issuing 

from the tank (neglect the kinetic energy at section 1)? 


Solution: 


The energy equation for section 1 and 2 


2 2 
jy iecE lg. ! =z,+P24 we 
pg 2g pg 2g 


Use Section 2 as the datum, so 
2 


h+0+0=0+0+—2 
2g 


Therefore 


u, =./2gh 
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2) Energy Loss and Gain 


In practice, some energy is ‘lost' through friction (hf), and external energy may be added by means of a 


pump or extracted by a turbine (E). The energy equation will be 


Py 


ue 
ztt+t+E=z,+ 
PS 


& 


3) Venturi Meter 


is an instrument which may be used to measure discharge in pipelines. By measuring the difference in 


pressure, discharge may be made as (assume no energy loss) 


2 
eon eae ; 
pg 2g 
i.e. 
a ae 
pg 2g pg 2g 


From the continuity equation 


Q= Au, = Au, 


A, 


V(A/ A)? -1 


pg 2g 


Section 1 
Section 2 


The actual discharge will be slightly less than this due to the energy losses between Sections 1 and 2. 


A coefficient of discharge is introduced to take account of these energy loses. 


Q cual =C gQ (Cd values can be found in British Standard 1042 or other hydraulic handbooks) 


Orifice plate 


is similar to Venturi meter with compact size (hence more economical) 
but less accurate and more energy loss. 
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Questions 3 
Energy Equation 
1. The device Pitot tube shown below is used to determine the velocity of liquid at point 1. It is a tube 
with its lower end directed upstream and its other leg vertical and open to the atmosphere. The impact 


of the liquid against Opening 1 forces liquid to rise in the vertical leg to the height Az above the free 
surface. Determine the velocity at 1. 


(Answer: u = J2gAz) 
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2. What should the water level h be for the free jet just to clear the wall? 


50 cm 
Thin wall 


(Answer: 0.156 m) 


3. If the velocity at Point A is 18m/s, what is the pressure at Point B if we neglect friction? 


A 
D=75mm 
15m | 
D=200mm 
as 
0.5m 


——— Water 


(Answer: 308KN/m’) 


Solutions 3 


Energy Equation 


1. Kinetic energy at 1 is converted into Az increase of the fluid at 2 


2b Th 
Az 
2 
eee = 
2g a 
hence 
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u,=~2gh 


Fall distance gt~/2 = 0.40 
So t=0.8944/,/g 


Horizontal distance = ut = (f2 gh)(0.8944 / V2) = 0.50 


h=0.156m 


3. The energy equation between A-B 


V2 y 
zpt et eaz, +e 
pg 2g pe 2g 
i.e. 
es 13° 
0422 4-2 =15+0.54+—4+— 
pg 2g pg 2g 


p,/p=—0.5V,; +3141 


The energy equation between A-C, 


V- V- 
2g pg 2g 
i.e. 
0 Vv 0 18° 
0-4 —— = 154+ —— 
pg 2 pg 2g 
V. = 24.86m/s 
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4. Momentum Equation 


4.1 The Principle 


From Newton’s 2" law 


F=ma 
where a= 2— Vi 
At 
m(V,-V,) =m 
So ee 


where PQV is the momentum, i.e. Momentum = Mass x Velocity 


It is important to note that momentum and velocity are both vectors. 


atter what you want out of your future 

per, an employer with a broad range of 
operations in a load of countries will always 
be the ticket. Working within the Volvo Group 
means more than 100,000 friends and 
BSiieagues in more,than 185 countries all 
over the world. We offer graduates great 


career opportunities.- check out the Career 
section at our web site www.volvogroup.com. 
le look forward to getting to know you! 


VOLVO 


AB VoWo (pub!) 
Wr YONOGTOUD. COM 


Vowo Teucxs | Rexauer Teucns | Macx Trucxs | Vouo Buses | Vowo Coxsteucrion Equtpwext | Vowo Pena | VowoAeno | Vowo IT 
Vowo Financut Serwnces | Vowo 3P 1! Vowo Powerman | VowoPaers | Vowo Tecumoocr | VowoLoasnes | Business Anca Asw 


Download free ebooks at bookboon.com 


32 


Concise Hydraulics Momentum Equation 


Therefore, momentum force is a vector: F = pQ(u, —U,) 


It is more convenient to use scalar forms by separating the momentum force M y 


into three basic components: 


E = PQ(Uu,, — Ux) 
f = PQ(Uyy — Uy x 
i = pQ(u,, =i) 


4.2 Applications 

Two steps: 

a) Calculating the total momentum force; b) find out where this force is sourced. 
1) A jet normal to a fixed plate 

Estimate F'R ( the force exerted on the fluid by the plate) 

Solution: 

a) In the horizontal direction, the total momentum force is 

—F = pQ(u,, —U,,) 

Here 

u 


=0,,u,.= u, Q=Au 


2x 
F = pAu’ 


b) Since there is only one force FR acting on the water jet, the total momentum force equals to FR 


F, = F = pAu* 
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2) Force exerted by a nozzle 


Calculate the force FR required to hold a nozzle to the firehose for a discharge of 5 litre/second if the 


nozzle has an inlet diameter of 75 mm and an outlet diameter of 25 mm. 


Solution: 


a) The total momentum force can be derived from 


Fp <+——_ 
P 
—> —= 
uy u2 
2 
1 
Control volume (Grey area) 
F = pQUu, —u,) 
2 
A= —— ) 4.42102 m? 
2 
= P0025) _ void ta? 


u, =Q/A =1.131m/s 
u, = Q/A, =10.18m/s 


So the total momentum force on the fluid is 


F = pQ(u, —u,) = 1000x 0.005 x (10.18 —1.131) = 45.25N 
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b) 
Force balance 
F=pA-F, 


From the energy equation 
u, ob P, = u, a P» 
2g pg 29 pg 


Therefore 


2 2 
p = PO A) 51.2kN fm 


1 


Substitute into 
PepA =f, 


F, = p,A,—F =51.2x4.42— 45.25 =181N 


What do you see? 


A skateboarder? 


We see an All Options Junior Trader leaving 
the office after a successful day. 
He has the freedom to make his own decisions. 
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3) Two dimensional flow 


Calculate the magnitude and direction of the force exerted by the pipe bend if the diameter is 600 mm, 
the discharge is 0.3m3/s and the pressure head at both end is 30m. 


Solution: 


2 
a) A= as) = 0.283m? 


u= Q/A=1.06m/s 
So the total momentum force on the water 


F, = pQu,, ~U,,) = PQ(0-u) = -318N 
F, = pQ(u, -u4,)= pQu - 0) = 318N 


b) The total momentum force is a resultant force of pressure force and reaction force from the pipe. 


F,=F,,+Ap and F, = F,,— Ap 
F,, = F, — Ap = —318 — 83202 = —83.52kN 
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Fe, 


y 


= F, + Ap = 318+ 83202 = 83.53kN 


So 


F.= VF. + Fa, =118.1kN 


and@ =tan ‘(Fy / Fa.) = 135° 


Questions 4 
Momentum Equation 


1. A small ingot and platform rest on a steady water jet. If the total weight supported is 825N, what is 
the jet velocity? Neglect energy losses. 


mm 


Water Jet 
D=60mm 


(Answer: 17.1m/s) 


2. Calculate the magnitude and direction of the force exerted by the T-junction if the water discharges 
are Q,=0.3m"/s, Q.=0.15m*/s and Q;=0.15m*/s, the diameters are D,;=450mm, D.=300mm and 
D;=200mm, and the upstream pressure p,=500kN/m’, Neglect energy losses. 
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3. A vertical jet of water leaves a nozzle at a speed of 10m/s and a diameter of 20mm. It suspends a 
plate having a mass of 1.5kg. What is the vertical distance h? 


(Answer:3.98m ) 
Solutions 4 
Momentum Equation 
1, 
70.0607 
A=—— =2.826m’ 


From the momentum equation 


F = pQu, —u,) =—pu,Au, =—2.826u; 


So 

2.826u, = 825 
u, =17.1m/s 
2 


A,= 0.159, A, = 0.07065, A, = 0.0314 
u, = Q/ A, =1.886m/s, u, =2.123m/s, u,=4.775m/s 


Pressure 
[a ee ae 
p9 2g pg 2g pg 2g 


2442 
P, = p+ Aa = 499.SKN / m? 


Download free ebooks at bookboon.com 


39 


Concise Hydraulics Momentum Equation 


2 9 
Ps=D, +o) = 490.4KN / m? 


The total external force 
F. = P(Quu,, + Qus, — Qu,,) = —565.8N 
EB, = P(Q,U, - Quus, - Quy) =—-397.8N 


Reaction force 


Fo= PAT FA PA Pit Fy 


F., = FE — p,A, =—565.8 — 500 x 159 = —80.1kN, 
y = F, — pA, + pA, = -397.8 — 490.4 x 31.4 + 499.4 x 70.65 = 19.48kN 


Hence 


Fa = [Fe + Fg, =82.43kN 


@ =tan”' (Fry / Fe, )= 166.3° 


Be 

The required speed to suspend the plate 
Q=VA=10x 0.02* 7/4 = 0.00314m’ /s 

F. = pQ(u,, —u,,) =—1000 x 0.00314V = 1.5g 


So 
15g 


1000 x 0.00314 


The energy equation 


2 2 
Pry 4M a Prag gt 
pg 2g pg 2g 
2 2 2 2 
heap We = 3.98m 
2g 2g 2g 
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5. Pipe Flow 


5.1 Introduction 


Pipes are widely used in engineering to deliver fluid from one place to another. There are two 
flow types. 


Laminar : flow in discrete layers with no mixing Laminar 


Turbulent: flow with eddying or mixing action. 


; Sh) Turbulent 


Reynolds Number is used to predict the laminar/turbulent flows 


ey 
a 


Re 
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The Reynolds number represents a ratio of forces 

inertia force 
Re = ———— 
viscous force 


Re < 2000, laminar, 
Re>4000, turbulent, 
2000<Re<4000, transition region 


Velocity distribution in laminar flow 


1s 
U=Unax ~ R2)? 


where U.,,,, = 2U 


Velocity distribution in turbulent flow 


( 


1/7 
u= Vinax | 1 — 4 , where u,,,, = 1.23u 


max 


Mean velocity 


Vag = 


The energy and momentum equations may then be rewritten in terms of the mean velocity 
av? 
Oe and F = pQftV, -V,) 


In most pipe flow, ~=1.06 and #=1.02 and both may safely be ignored in practice. 


5.2 Energy losses in pipe flow 


Darcy-Weisbach equation 


1 2 


where hf is the energy head loss, L is the pipe length, D is the diameter and V is the mean velocity. 
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The friction factor / is dimensionless and can be derived by experiments or from the following chart 
or formulas. 


For laminar flow 


For turbulent flow (Colebrook-White 
formula) 


(_k 351) 


Ss 


1 
a 54 + 
a °F3-7D " Reva/ 


Ks is the roughness factor 


Typical Ks values (in mm) 


The Moody Chart 
brass, copper, glass 0.003 
plastic 0.03 
galvanised iron 0.15 
Example 


Oil flows through a 25 mm diameter pipe with a mean velocity of 0.3 m/s. Given that 1=4.8x107 
kg/ms and p=800kg/m”, calculate (a) the pressure drop in a 45 m length and the maximum velocity, 
and (b) the velocity 5mm from the pipe wall. 


Solution: 


1) Check if the flow is laminar or turbulent 
Re = pDV/ 1 =800x0.025x 0.3/(4.8x10~) =125 < 2000 


So it is laminar flow 
2) Pressure drop 


A= = = 64/125=0.512 


Re 
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AV? 0.512x45x0.3" 
h, ee Aon 
D2g  0.025x2g 


or convert it into kN/m? 


pgh, =33.2KN /m* 


3) Velocity 


Unax = 2V =0.6m/s 


u = 0.384m/s 
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5.3 Local losses (or Minor losses) 


Local losses occur at pipe bends, junctions and valves, etc. 


Vv? 
Lao 


h, =K 
L 2g 


where /h, is the local head loss and K; is a loss factor for a particular fitting (usually derived from 


experiments). V is the mean flow velocity before or after the pipe fitting (use the higher value if they 
are different). 


The energy equation for a complete pipe with fittings will be 


2 2 
z tebe tai tye h, +h, 
2 Ps& & 


5.4 Grade Line 


Hydraulic grade lines: A line that represents the potential energy and pressure energy in the fluid. 


ees 


pg 


When a hydraulic grade line is below the pipeline, the pressure in the pipe is less than the atmospheric 
pressure (negative pressure). 


Energy grade lines: A line that represents the total energy in the fluid. 
P 2 


i a ry 
pg 2g 


It is always vertically above the hydraulic grade line by a distance of ve /(2g). 
5.5 Combination of pipes 
1) Pipes in series 


Two or more pipes of different sizes or roughnesses are so connected that fluid flows through one pipe 
and then through the others in turn. 
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H= > local losst+hy + hy, 


Q, =Q, 


2) Pipes in parallel 
The flow is divided among the pipes and then is joined again. 
Q= Q, + Q, + Q; 


Energy loss in pipe 1 = Energy loss in pipe 2 =... 


3) Networks of Pipes 


Problems on networks of pipes in general are complicated and require trial solutions in which the 
elementary circuits are balanced in turn until all conditions for the flow are satisfied. 
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Binal 
| L 


5.6 Energy Loss in Non-circular Pipes 


The Darcy-Weisbach equation for circular pipes can also be used by replacing diameter D with 


hydraulic radius R. 


_ Area 
Perimeter 


It can be applied to square, oval, triangular, and similar types of sections. 
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Questions 5 
Pipe Flow 
1. In the laminar flow of a fluid in a circular pipe, the velocity profile is exactly a true parabola. The 


discharge is then represented by the volume of a paraboloid. Prove that for this case the ratio of the 
mean velocity to the maximum velocity is 0.5. (Note: y= Unasl —(r/ R) | ) 


2. An oil (p=868.5 kg/m? and =0.0814 kg/ms) is to flow through a 300m level concrete pipe. What 
size pipe will carry 0.0142m3/s with a pressure drop due to friction of 23.94kPa? 


(Answer: 155 mm) 


3. Water flows at a rate of 0.04m’/s in a 0.12-m-diameter horizontal pipe that contains a sudden 
contraction to a 0.06-m-diameter pipe (K,=0.40). Determine the pressure drop across the contraction 
section. How much of this pressure difference is due to energy losses and how much is due to kinetic 
energy change? 


(Answer: 133.9 kPa , 40.0 kPa, 93.9kPa@) 


4. If H = 11m, find the discharge through pipes 1, 2, and 3 (Neglect local losses). 


Pipe 1 70m D=50mm, A=0.114 
Pipe 2 80m D=120mm, ,=0.088 
Pipe 3 110m D=100mm, ,=0.114 


(Answer: 0.0009, 0.00856, 0.00946 m3/s) 


Solutions 5 
Pipe Flow 


1. 


Q=[udA ={ ya} - = mur 2 tal = 


Thus V/u,,, = 0.5 
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D 
ras =0.0142, D=-+40.0181/V =0.135m 


Trial 1: V=1 m/s, from Q = AV = 


Re = —— = ———"—— = 1440 (laminar) 


4 | Res 
Bees =0.0444 and h, = js te high, 
R , D2 


€ § 


LW 
Sub hf with 2.81 m to estimate V. h,= Crs = 2.0m 
& 


V=0.75Im/s 


Trial 2 V=0.751 m/s soD=-+0.0181/V =0.155m 


gt Se EE? 194) anal) 


64 LV 
A=— =0.0515 so h, = aie =2.87m close to 2.81 m 
R D2 


€ § 
Hence a pipe diameter of 0.155m would be required. 


3. From energy equation 
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Pi- Pp _Vi-Vi , ag Vy _1415°-354" | 14.15" 


pg 2g "2g 2g 2g 
= 9.57 +4.08 = 1365m 


P, — P2 =93.9 + 40.0 = 133.9kPa 


This represents a 40.0 kPa drop from energy losses and a 93.9kPa drop due to an increase in kinetic 
energy. 


4. 

0; QO, QO, Q, Q, a, 

==! = —=_] =509.60,, V, == = — = 88.50,, V, == =—5 = 1274 
' A, 70.0257 Qu ts A, 70.067 22s A, 70.05” e; 
hint 
i.e. 
ee Ve L, Vy 


600T 'A'g SweIsAS WIy] JEW] @ 


It’s only an 
opportunity if 


you act on it 


IKEA.SE/STUDENT 
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So 

chy ee 70 (509. 6Q,) ~ 0.088 80 (88. (88.5Q,)* 
0.05 22 0.12 2g 

Q, =0.10590, 


Q0,=2,+92,=1.1059, 
2 2 
jig 
D, 2g D, 2¢ 


80 (88. (88.5Q,)° 110 (127.4Q,)” 


11 = 0.088—— +0.114— . 
0.12 0.1 2g 
oe iO gpg cp ai Oey 
012 2g 0.1 2g 


Q, =0.00856m'/s, Q, =0.105Q, = 0.0009m*/s, QO, = 1.1059, 


= 0,00946m’/ s 


Download free ebooks at bookboon.com 


Dl 


Concise Hydraulics Physical Modelling 


6. Physical Modelling 


6.1 Background 


Hydraulic modelling includes physical and mathematical models. Initially, due to a lack of computing 
power, physical models dominated the hydraulic modelling field. Froude (1870) established the famous 
Froude similarity when he was testing a ship model. In 1885, Reynolds built a river model for Mersey 
based on the Froude similarity. The Buckingham’s Pi theorem (by J. Buckingham in 1914) extended 
the similarity principle to a broader sense. Since then, physical models have been widely used in 
hydraulics, including complex 3D river models with sediments. 


With the advent of computers, mathematical modelling has been gaining popularity. It is possible to 
solve complicated differential equations using numerical techniques due to modern computer’s large 
memory and CPU speed. Nowadays, mathematical models are highly accurate and efficient to tackle 
1D river problems. They are cheaper to build and fast to run in comparison with physical models. In 
2D modelling, mathematical models are gaining ground from physical models. However, for 3D 
modelling, physical models still have superior performance in many complicated problems and it 
would be a long time before mathematical models are able to totally replace them. 


6.2 Dimensional Analysis 


Dimensional Analysis is a branch of applied mathematics for investigating the form of a physical 
relationship. Dimensional Analysis not only helps us to design experiments but also allows us to gain 
an insight into how all the factors depend on each other even before we start experimenting. 


Example: 


A thin rectangular plate having a width w and a height h is located so that it is normal to a moving 
stream of fluid. Assume the drag D, that the fluid exerts on the plate is a function of w, h, the fluid 
viscosity t1 and density p, and the velocity V of the fluid approaching the plate. Perform a 
dimensional analysis of this problem. 


WwW 
Basic steps: 
LH, p 
1. Decide variables to be used SS 
V h 


D= f(w,hu pV) 
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2. Express each of the variables in terms of basic dimensions 


mass M, length L and time T (or force F, L, T) 


w=L,h=L,u“=ML'T" (= dimensionally equal to) 


p=ML*, V=LT", D=MLT? 


3. Determine the required number of /Tterms (free variables) 


Buckingham [T theorem: the number of /Tterms (free variable) is equal to k-r, where k is the total 


number of variables in the problem and r is the number of reference dimensions (up to 3). 


k-r = 6-3=3 so we have three free variables 


4. Select dependent variables, 


Select w,V and p as dependent variables ( Avoid choosing both w and h since they have the same 


dimension). The rest are free variables (D,u,h). 


<< 
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5. Form a [Tterm by grouping one of the free variables with all dependent variables, 


Starting with the free variable, D, the first I] term can be formed by combining D with the dependent 
variables such that 


ay7b c 
Il, =Dw Vp 
Estimate a,b,c so that the combination is dimensionless 


(MLT ’)(L)*(LT')?(ML*)° = M°DT® 


So 

1+c =0 (for M) 
1+a+b-3c=0 (for L) 

-2-b =0 (for T) 
Therefore 


a= -2, b =-2, c= -1, the IJ term then becomes 


6. Repeat step 5 for each of the remaining free variables, 
Il, = hw°V"', and a=-1, b=0, c=0 

So II, =h/w 

And 


TI, = uw'V" p*, and a=-1, b=-1, c=-1 


So ,=—— 


wVp 
7. Check all the resulting /7terms to make sure they are dimensionless, 


hk 
Il, =—=—=L’, the same for the other two 
w L 
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8. Express the final form as a relationship among the /7terms 


Dd rE _) 
w Vp w' wVp 
Now, there are only three variables instead of six, which are much easier for experimental data analysis. 


Dependent variables are also called repeating variables, and free variables are called nonrepeating ones. 
6.3 Analysis of Experimental Data 


Dimensional analysis can greatly reduce the number of experiments needed and improve the data 
processing efficiency. 


1) One [Iterm 

Example: Assume that the drag D, acting on a spherical particle that falls very slowly through a 
viscous fluid is a function of the particle diameter d, the particle velocity V, and the fluid viscosity u.. 
Determine how the drag depends on the particle velocity. 

Solution: 


D= f(d,V, 1) 


by dimensional analysis, only one /7, II, = Constant 


D 
=; 
iWvd 
So 
D 
——=¢ 
LVd 


C was found by experiment C=3z, so 
D=3z2uWVd (Stoke law) 


Therefore, only one experiment is needed to solve this problem. 
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2) Two II terms 

= g(t.) 

For problems involving two I terms, results of an experiment can be conveniently presented in a 
simple graph. The relationship between the two terms can be solved with curve fitting techniques (e.g., 
linear or nonlinear regression). 

3) Three IJ terms or more 

= g(11,,11,) 

It is still possible represent three variable problems in 3D graphics. Mathematically, surface fitting is 


needed to solve the relationship among them. If more than three terms are involved, more complicated 
data mining techniques are needed (such as Artificial Neural Networks, Support Vector Machines, etc). 
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6.4 Model and Similarity 


A model is a representation of a physical system that may be used to predict the behaviour of the 
system of interest. 


As we know, any given problem can be described in terms of a set of I] terms as 


1 SO eu, Kt ) 
For a model to represent such a system 
Tin = (TI, jl ies aig) 


Therefore 


if  ,,=H,,H,,,=0,,...I1,, =, then I, =I1,, 


fo ae nm n 


This is the prediction equation and indicates that the measured value of I1,,, obtained with the model 


will be equal to the corresponding II, for the prototype as long as other pi terms are equal. 


In many hydraulic engineering problems, models are with free surfaces (rivers, estuaries etc). Force F 
is linked with length, velocity, viscosity, and gravity acceleration, hence f(F,1,V,p,u,g)=0 . 


From which, dimensional analysis gives: 


F VI V° F 
72 =f alae or eo = f (Re, Fr), 
pv'l u gl pv'l 


l 
where Reynolds number Re = ee , Froude number Fr = = (to be introduced in later chapters). 


“ Jal 


In practice, it is usually impossible to satisfy all the ,] terms, so it is important to identify the 
dominant terms. The Froude number is very common in river modelling, which represents 
gravitational forces in dominance. Hence, the Froude numbers for the model and the real-world 
structure (prototype) should be the same. 
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Froude Number Fr, = Fr, =V,/./gl, =V,,/<(gl, » hence V,,/V, =,/L,,/L, = JA, 


Similarities 

Geometric Kinematic (velocity) Dynamic (force) 

cP L_/T F OM ic 

fea, Vent Ett Palast aga A) ait 
i Vi ytd, we B. Mid. pb, Ae A, 


Questions 6 
1. Water flows over a dam as illustrated below. Assume the flow rate q (per unit length along the dam, 
hence with unit of m’/s) depends on the head H, width b, acceleration of gravity g, fluid density p, and 
fluid viscosity 1. Develop a suitable set of dimensionless parameters for this problem using b, g, and p 
as repeating variables. 


q AH 
(Answer: 73/2 1/2 = (ana )) 
b’’g b b''g"p 


Hl aq—> 
1 = 


ae ad 


2. A 1:50 model of a boat has a wave resistance of 0.02N when operating in water at 1.0m/s. Find the 
corresponding prototype wave resistance. Find also the horsepower requirement for the prototype. 
What velocity does this test represent in the prototype? (use Froude criterion, 1Kilowatt = 1.34 
Horsepower) 
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Solutions 6 
Physical Modelling 


1. 
q= f (Lb, 9, p, 4) 
_f ae oy _L .M 
ea —o Lig poP > a3? ws ML"T* 


Six variables, only three reference dimensions, 
so k-r = 6-3 =3, three IT terms 


Repeating variables: b,g and p 


HL, = qb‘g’ p* 2s Cry Lr) ae = [°T°M° 
2+a+b-—3c=0 

—1-2b=0 

c=0 
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Therefore c=0; b=-1/2; a=-3/2 


So, Th = pan giz 


Il, =Hb*g’o° =L(L)°C.T-y (Mey = PT’ 


1+a+b-—3c=0 
—2b=0 
c=0 
Therefore, 
c=0; b=0; a=-1 
n= 2 

b 


N,=(6'¢ o = MET ay ary (Mey arr Mt 
—l+a+b-3c=0 

—1-2b=0 

1+c=0 


Hence c=-1; b=-1/2; a=-3/2 


Lu 
I], = bg! 5 
q es 

Therefore bi2gi? = 9( 5° be gZp 
2. 
The geometric ratio is 

EL 1 
A, = =—=0.02 

L 50 


P 


From the Froude number 


Froude number Fr, = Fr, =V,/,/gl, =Vn/ 9m » 
hence Ay =V,,/V, =4/L,,/L, = 4, =V0.02 = 0.1414 


The velocity of 1.0m/s in the model represents the prototype velocity 
V, =V, {0.1414 =1/0.1414 =7.07m/s 
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Dynamic (force) similarity can be derived as 


F =mass x acceleration 


F, M LT digi, A 002 
m mdm _ Pm ann di 4 = ue (The fluid density is the same in both cases) 
F, M,a, p,L, 1,1, A; A; 


Pp 


The time ratio can be derived from Kinematic (velocity) similarity 


V, L,/T, A; 


m — = =A, 
4 ee A, 
Hence 
fees eee =0.1414 
Ay 0.1414 
So 


2 2 
oe — Bi — 0.02 = 125000 0.02 = 2500 


The force ratio is 


m 


F, 125000 


F 


From Power =V x F 


P, = FV, = 2500 x 7.07(watt) = 17.68(kilowatt) = 23.67hp 
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7. Open Channel Flow 


7.1 What is "Open Channel Flow" 


Open channel flow is flow of liquid in a conduit in which the upper surface of the liquid is in contact 
with the atmosphere. The flows in the following cases are all open channel flow (include the one in the 


There are many types of open channel flow and a common classification system is as follows. 


Classifications 


circular pipe). 


Steady Unsteady 
Flow Flow 


Uniform Varied Flow 
Flow 
Gradually Varied Rapidly 
Flow Varied Flow 


The uniform flow refers to flow whose water depth, width, flow area and velocity do not change with 
distance. If any of these factors change, the flow will be called varied flow (or nonuniform flow). 


=o 
D, 
| D 


2 


z Ey 
Diep, Ditty 
Cy Y 
CZ Varied flow 7, 


Uniform flow 


Constant depth 
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The steady flow refers to flow factors and variables that do not change with time, otherwise the flow 
will be called unsteady flow. 
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7.2 Channel Geometric Properties 


The measurable geometric properties of channels are: 
Depth (y), Area (A), Wetted perimeter (P), Surface width (B) 


They are aggregated into two useful factors: 
Hydraulic radius (R=A/P), Hydraulic mean depth (Dm=A/B) 


-— 8 
4 |} —s—4, 
».« 
kK— b—4 4 <— p—_> # 
Rectangle Trapezoid Circular 
i by (b + xy)y ; (¢—sing)D° 
p b+2y b+2yv1+ x? 50D 
B b b+2xy (sin) p 
R by MOP y _ i sing) 
b+2y b+ 2yv1+ x? 4 = ~p ye 
: ’ (b+ xy)y 1( g=sing) 
m b +2xy 8 \ sin(g/ 2) 


Irregular channels 
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A= “(h-y)dx P=| "dx +dy’ =[* 1+ y"dx 


In practice, the cross section is divided into a finite number of sections. A and P are the sum of those 
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7.3 Calculation of Hydraulic Radius and Hydraulic Mean Depth 


Example 1: Rectangular channel 


1m 


4 


|——3m —| 


From A=3m’,B=3m, P =1+1++3=5m 
So R=3/5=0.6m, D, =3/3=1m 


y 
x _/* 
45 


fe 3m 4} Z. 


Example 2: Trapezoidal Channel 


From A= 4m’, B=5m, P =5.83 
So R=4/5.83=0.686m, D, =4/5=0.8m 


Example 3: Circular Channel 


eae 4 


From B= 2 =1.414m, P =3.14/2=1.57m, A=3.14/4-1/2 =0.285m* 
So R=0.285/1.57 =0.182m, D,, =0.285/1.414 =0.202m 
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Questions 7 
Open Channel Flow 


1. Compute the hydraulic radius and hydraulic mean depth for a trapezoidal channel. 


(Answer: 0.91m, 1.2m) 


2. Compute the hydraulic radius and hydraulic mean depth for a smooth concrete-lined channel. 


(Answer 0.735m, 1.0 m) 


3. Derive equations for hydraulic radius and hydraulic mean depth for a parabolic channel. If B=6 m 
and h = 3m for the channel, calculate its hydraulic radius and hydraulic mean depth. 
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P (B/2, h) 


0 
Function for a parabolic curve y=kx? 


Note: 


[ Vax’ + cdx = 5 Vax” +¢€ +575 Inova +Vax’ +c) 
a 


(Answer 1.35m. 2.0m) 
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Solutions 7 
Open Channel Flow 
oo) L5 
_ 6+3)x15 _. 2 
2 
P=3+42V1+1.5" =6.61m 
Hydraulic radius R= A/P=6/6.61= 0.91 m 
B=5 m 
Hydraulic mean depth Dm=A/B = 6/5= 1.2 m 


1. A 


2; 


A=3x15+ 2" = 6m" 


P=34+154+054++714+3? =816m 
B=6m 


Hydraulic radius R=6/8.16=0.735m 
Hydraulic mean depth D,,=6/6=10m 
3. 


Assume the parabolic curve function is y=kx2 


bg 


(B/2, h) 


Replace (x,y) with (B/2,h) 
So k= 4h/B2 


The following steps can be done either manually, or by Matlab (or Maple) 
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Method 1 (manual) 
B 2 
Area A= J2,(h — y)dx = aul 
2 


B B 
P =J2pylty dx = 2)2y14 y" dx 
“2 


B 
= 2[3 {1+ (2kx)’ dx 


Since ax” +cdx = ax” +c +— = In(xva +-Vax° +c 
Nat scars or eee ine ela 40) 


2Va 


= 
2 
and let s=(2k)2 


B Ba 
P = 2jav1+ sx’ dx = 25 Asx" +1 + spines +¥sx’ +1) |?) 


Bila. Bvs {BY 
5 (5) +147 as 5) +1) 


substitute s with (2k)2 


2 ew'(2) ) 1+ ine + yen) +1) 


- 2 aay #19 5, in(kB + (kBY +1) 
SRC) 1+ Sink + (0 FD 


let t=kB = 4h/B, then 


P=(B/2)W1i+e +1/thn(t+Vi+e)] 


Since 
R=A/P 
so equations for hydraulic radius and hydraulic mean depth are: 


4h 


Ah 
R ae where t = — 
3WW1 +0? +1/tin(t+vi+e)] B 
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and 


4h 4x3 
t=— =— = 2 
B 6 
Fe 
3[W14t? +1/tIn(t+-V¥1+t?)] v1 4 22 +tin(2 +142") 
2 
= 4 = —_~$ «1.35 
[/5+ 5 in +75)] [2.236 +50.444)] 
3 3 
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Method 2: Matlab (it actually uses Maple as its symbolic engine hence you can also use Maple 


directly) 


Write the following lines into openchannel.m, then run it: 


%-- Matlab code for Open Channel Q3 


syms x kh B 

%-- curve equation 
k=4*h/B/2; 
yHk*x 2; 


%- area A 
A=int(h-y, -B/2, B/2); 


% wetted perimeter 


P=int((1+(diff(y))/2)"0.5,-B/2,B/2); 


%-- hydraulic factor 
R=A/P; 
R=simple(R) 
Dm=A/B 


%--- substitute values B=6, h=3 
B=6;h=3; 

eval(R) 

eval(Dm) 


Results: 
R — 


32*h\2*B/(24*(BA2+ 16*h\2)\(1/2)*h+3*log((h\2*64\(1/2)+2*((BA2+ 16*h\2)/BA2)\(1/2)*BN3*(h\2/ 


BM)N1/2))/B13/(W\2/BM4)\(1/2)) *B2-3*log((- 


h\2*64\1/2)+2*((BA2+ 16*h\2)/B/2)\(1/2)*BA3*(h\2/BN4)\1/2))/BN3/(H\2/BM4)\(1/2))*B2) 


Dm = 
2/3*h 
R= 
1.3523 
Dm = 
2 
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8. Uniform Flow 


8.1 Introduction 

Definition: Water depth, width, flow area and velocity do not change with distance. 

An important term in hydraulics is flow rate or discharge. It is defined as the volume of water that 
passes a particular reference point in a unit of time. Its common units are m’/s or cumecs (other units: 
cfs , mgd, etc). 

To understand how large 1m”/s is, here is an example. 

Example: 

If the daily water consumption in Bristol is 300 liters/person, how many people would be provided 
with water by 1 m°/s flow? 

Solution: 

1000 x 24 x 3600(I / day)/300 = 288, 000 

So it could supply water to 288,000 people in Bristol. 


8.2 Laminar or Turbulent Flow 


As in a pipe flow, the Reynolds Number (Re) can be used to identify laminar and turbulent flow state 
for open channels. 


RV 
ee al 
L 
where : 
R hydraulic radius (m) p density (kg/m), 1000 kg/m? 
Vv mean velocity (m/s) wm viscosity (kg/m s), about 1.14 x 10° Ns/m? 
Laminar Flow Re < 500 
Turbulent Flow Re > 1000 


For practical open channel flows in civil engineering, usually | Re >> 1000, i.e., turbulent flow (try 
Re for R=1m and V=1m/s). 
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8.3 Energy Loss Equations 


What do you see? 


A skateboarder? 


We see an All Options Junior Trader leaving 
the office after a successful day. 
He has the freedom to make his own decisions. 
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Rearrange it into 


So 
V=CyRS, where S,=h,/L,and C=./2g/A 


This is the famous Chézy Equation and C = the Chézy coefficient. C is not constant and depends on the 
Reynolds Number and boundary roughness. 


Another derivation with the coefficient by setting C= R’’°/n results in the Manning Equation. 
Manning’s Equation is the most widely used formula in engineering and all the examples and questions 
in this book are based on this formula. 


Manning’s Equation 


, _ R23 Sy 
n 
where n= Manning’s coefficient 


Since Q=AV, another form of the Manning Equation is 


ATS, 


Q _ npee 


Manning's n can be estimated by measuring Q, A, P, and Sg for a chosen channel reach & cross section. 


In practice, many researchers have carried out a large number of experiments to estimate Manning’s n 
and they can be found in most hydraulics handbooks. 


Typical values of Manning's n 


Channel type Surface n 

River earth 0.02 - 0.025 
gravel 0.03. - 0.04 

Unlined canals earth 0.018 - 0.025 
rock 0.025 - 0.045 

Lined canals concrete 0.012 - 0.017 


Note: More values can be found in Chow’s Open Channel Hydraulics (1959). 
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8.4 Computation of Uniform Flow 

There are two common problems to solve: 

1) Determine the discharge given the depth 

Example: 

The normal depth of flow in a trapezoidal concrete lined channel is 2m. The channel has a base width 
of 5 m and side slope of 1 (vertical):2 (Horizontal). Manning's n is 0.015 and the bed slope Sg is 0.001. 
Determine the discharge Q and the mean velocity V. 

Solution: 

Known: b=5m, y=2m, x=2, n=0.015, So=0.001 


Unknown: Q, V 


A=(b+xy)y =(5+2x2)x2=18m* 


P=b+2yV1+x? =54+2x2V5 =13.94m 
R=A/P=1.29m 


Compute V first, then Q 


RPS; 1.99°* 001 


n 0.015 


V =2.5m/s 


Q = AV =18x2.5=45m?/s 
2) Determine the depth given the discharge 


Example: 
If the discharge in the channel given in the last example is 30 m’/s, find the normal depth of the flow. 


Solution: 
Known: b=5m, Q=30 m’/s, x = 2, n= 0.015 and So=0.001 
Unknown: y 


A= (btxy)jy=(5+2y) y 


P=b+2yVl4¢x =54+2v5y 
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AX3YS,  [(5+2y)y]_ x-V¥0.001 
nP? 0.015 x [5+ 2V5y] 


5/3 
di 14.23 = Ot 2y)yT 
(5+ W5y)"’ 


Trial-and-error method (try to find a y value which makes the Right-Hand-Side close to 14.23) 


y RHS 
1 5.72 
1.5 12.20 
1.6 13.80 
1.7 15.52 
1.65 14.65 


Hence the answer is y =1.65 m 
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Note: 
Slope description options (different descriptions for the same slope): 


1000 


Channel bed slope So= tan 0 
a) 0.001; b) 0.1 percent; c) 1:1000; d) 1 (V):1000 (H); e) 1m drop over 1km length 


Channel side slope 
a) 1:2 ( as the channel bed slope); b) 1 (V): 2(H) ; c) x=2  (X is the inverse of slope) 


Questions 8 

Uniform Flow 
1.Water is in steady uniform flow through a finished-concrete channel (Manning's n=0.012) as shown 
below. If the bed slope is So= 0.0016, what is the discharge? 


(Answer: 102 m*/s) 


0.9m 


3.1m 


2. A trapezoidal channel has Manning’s n=0.020 and bed slope So=0.0004, which is made in the shape 
of a half-hexagon for maximum efficiency. For what length of L will the channel carry 7.1m°/s of 
water at its full capacity? 


60° 


(Answer: 2.33 m) 


3. A trapezoidal channel with a bottom width of 6m, a side slope of 1 (vertical):2 (horizontal), n=0.02, 
and a bed slope S,=0.0014, conveys a flow of 30m”/s. Estimate the normal depth of flow in this 
situation. 

(Answer: 1.60 m) 
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4. Given a circular culvert of 0.95 m in diameter with So=0.0016 and Manning's n=0.015, find the 
normal depth of flow for a discharge of 0.41 m’/s. 
(Answer: 0.511 m) 


Solutions 8 
Uniform Flow 
Is 
Area 
2 


mw-31 


A= +0.9 x 2x 31=20.7m’ 


Hydraulic radius 
A 20.7 


= — = ——— = 18m 
P 3174+2x09 


From the Manning Equation 


REVS, _» 9.7 18°" v0.0016 


Q=AV=A =102m’/s 
n 0.012 
2, 
ae 
L 
L 60° 
L 
A= Lsin60°(L+Lcos60°) = L’ — =1399i7 
P=3L 
R=A/ P=1299/3L =0.433L 
R*? S ; 0.667 | i 
Q=AV= aR vS. = 12991? oe)) ena 2 val 
n 0.02 
[287 = 955 
L=2.33m 
3. 


From the Manning Equation 
A= (b+ xy)y = (6+2y)y 


P= b+2yl+x = 6+ 2yv1 +2? =6+44.47y 
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pa bt2yy 
6+4.47y 


Ars 2y)yT’°>¥0.0014 
o MAYS. _[i6+2y wT" V00014 _ 44,9), 


nP? 0.02 x (6+ 4.47y) 
Simplifying 
Cyl aes 

(6+447y)> 

Try Left 

y=1 6.69 

y=1.5 14.05 

y=1.6 15.87 


So the normal depth is about 1.60 m 
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Known: D= 0.95 m, Q=0.41 m3/s 


n= 0.015, Sg=0.0016 


Unknown: y 


Using the following equations 


— Deng tl 4_ *y) 
Angle g@ = 2cos Ul 7 
1 ; 2 
Area = geen e)D 
1 
Wetted perimeter PS ; gD 


With Manning’s equation 


Le S 
gen 


np 
So 
0.41 - Le = sin g)0.95° /8]°°/0.0016 


0.015(0.95¢/2)°”° 


With n and So substituted 


Q 
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Pa Ps Pe 


0.046 


So $ =3.3, 
af =») 


From ¢@=2cos cars 


y = Dil -cos(@/2)]/2 = 0.95{1 — cos(3.3/2)]/2 =0.51m 


Method 2: Try different depths (slightly more calculations) 


-if 


2 
Guess a y and then use @ = 2 cos fe = =) to get the angle 


Soy =0.51 


Method 3: Matlab 
. 5/3 
Wilts 5o= (g—sin gy” 


g?'? 


Write a function in f.m: 


function y = f(x) 
y = (x-sin(x))\5/3)/Xx\(2/3)-3.65; 
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Channel Design 


9. Channel Design 


9.1 Channel Design 


Main purpose: to transport water between two points in a safe, cost-effective manner. 


Channel Types: Unlined , Lined(with concrete) , Grassed 


Permissible Velocities 


The maximum permissible velocity is the maximum velocity that will not cause erosion to the channel 


bed and walls. 


The minimum permissible velocity is the lowest velocity that will prevent both sedimentation and 


vegetative growth. In general, 0.6 to 0.9 m/s is used 
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Freeboard is the vertical distance from the top of a channel to the water surface at the design condition. 


Freeboard 


Water 


In design, 5% to 30% of the depth of flow is commonly used for freeboards. 


9.2 Compound Channel 


It is common to use compound channels in practical engineering projects to cope with different flows 
so that the sediments and vegetation in the channel could be controlled. However, the hydraulic 
computation for compound channels needs a special treatment as demonstrated by the following 
example. 


Example 1: 


E 


al 


3 m—>}+1 m>|<—3 m 
Case a: So=0.001, n=0.02, estimate the discharge in the channel 


Solution: 
The flow is contained in the main channel. 


A=by =1m’, P=b+2y =3m, R= A/P = 1/3 = 0.333m 


_R* JS, _ 0.333°°./0,001 


n 0.02 
Q = AV =1x0.76 = 0.76m?/s 


V =0.76m/s 


Case b: The same channel with a slightly more discharge (0.01m extra water) 
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Solution: 
y0.01m 


E 


al 
3 m—>}+1 m>|<—3 m 
The whole compound channel is treated as a whole cross section. 


A=1x1+7x0.01=1.07m? 
P=3+3+3+0.01+0.01=9.02m 
R= A/ P =1.07/9.02 =0.1186m 


R*,[S, _ 0.11867°./0.001 
n 0.02 
Q = AV =1.07x0.38 = 0.41m’"/s 


V= =0.38m/s 


The illogical result is due to the very different flow patterns between the main channel and the over 
bank flow. To solve this problem, it is necessary to treat different parts of the channel separately. 


Divided the channel into three parts: 


The main channel and two over-bank sections 


Areas : A1, A2, A3 
Wetted perimeter: Py(abc), Po(cdef), P3(fgh) 
Calculate V and Q separately and combine them together. 


Example 2 
During a large flood, the water level in the channel is given in the following figure. Manning's n is 


0.015 for the main channel and 0.035 for the flood banks. The bed slope S, is 0.001. Estimate the 
discharge for a maximum flood depth of 4 m. 
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i 14.5m > 15.0m it 14.5m id 


Solution: 


Split the section into subsections (1) in the middle, (2) and (3) as overbank parts and apply Manning’s 
equation to each one in turn and then sum the discharges. 


For Subsection 1 in the middle 


A, =15x1.5+(5+15)x 2.5/2 =47.5m’ 
P =5+(2V5x2.5) =16.18m 
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A’?/S, _ 47.5°°J0.001 


=—_V_9 = = 205.3m?/s 
Q np?" 0.015x16.187” 


Subsection 2 and Subsection 3 have the same dimensions, hence 


pie (Aa) x15 =18.38m’ 


P, =P, =10+(¥10 x1.5) =14.74m 


A/S, _ 18.38°°V0.001 _ 


3 
np? -0,035%14.74"" reer 


Q, =Q,; = 


Hence 
Q=Q1+Q2+Q3=243.7 m/s 
9.3 The Best Hydraulic Section 


A channel can have many shapes and dimensions. However, certain types of channels are better at 


conveying flows than others. 


Example: The following channels have the same amount of excavation (cross section area), as well as 
the same Manning’s n and bed slope. Which channel is the best hydraulically? 


(a) (b) 


Solution: 


Since P = 4, 4.24,8.5 


Ags K es 
np23 = pee (let KS n ) 


From Q = 


Then 
Q, =0.397K , Q, =0.382K , Q. =0.240K 


Hence a is the best and c is the worst. 
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A general pattern can be concluded from this example. 


ae A’? IS, 


From Manning's Equation = pie 


The channel having the least wetted perimeter for a given area has the maximum conveyance; such a 
section is known as the best hydraulic section. 


The derivation of the best hydraulic sections for some common shapes is illustrated here. 


Rectangular 
— = a 
y 
<a 
Ie b >| 
Let A = by = constant 
b= A/y 
So P=2y+b =2y+A/y 
Minimising P 
iP A 
dP = 2 —- —_ = 
dy y 
A b 
Then —>=2 or eo sO b=2y 
Ag Jy 


The best rectangular channel is the one with its depth being half of its width. 


i ty 
x 
kK b "| 


Trapezoidal 


A = (btxy)y 


P=b+2yV14+x* 


Now =a 


A 
Hence P=—~—xy+2yv1+x* 
J 
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OF 2 ao aay =0 
7. 


Oy 
oP 2x 
— + 


Solve both equations 


V3 


x =——=tan30° 
3 


and 
b=2y( lint 3} <2 22) 


Its is one-half of a hexagon 


k boo. 


6007 'A'g sueysAs yy] JeRU] @ 


R= It’s only an 
opportunity if 
you act on it 


IKEA.SE/STUDENT 
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Questions 9 
Channel Design 
1. Compute the discharge in a compound channel. The estimated Manning’s ns are 0.015 for the main 
channel and 0.035 for the over-bank sections . The bed slope of the channel S, is 0.001. 


(Answer: 97.09 m*/s) 


2. A rectangular, concrete lined channel (Manning’s n = 0.015) is to be constructed to carry flood water. 
The slope of the channel bed slope is 1 in 500. The design discharge is 10m’/s. (a) Calculate the 
proportions of the rectangular channel that will minimize excavation and result in the optimum 
hydraulic section. (b) If the cross sectional area of flow is kept the same as in part (a) but for safety 
reasons the depth of flow in the channel is limited to 1.00m, what will be the discharge now? 


(Answer: a) B=2y=2.89 m b) 9.55m’/s) 


3. Compare semicircular section with rectangular and trapezoidal sections to show that it is the most 
hydraulically efficient. 


4. Show that the best hydraulic isosceles triangular section is one-half of a square. 


Solutions 9 
Channel Design 
1. Split the section into subsections (1), (2) and (3) and apply Manning’s equation to each one in turn 
and then sum the discharges. 


For Subsection 1: B=3+2xy=3+2x2x2=11m, Ay=11x1+(3+11)x2/2 =25 m2 


Pp =3+(2yvi+x")=3+2x2Vi+4 =11.94m 
_ AS, 25°3.J0.001 _ 6.74 


ae Ne = 86.1m?/s 
Q nP?? ~—«0.015x11.9472_ 0.0782 


Subsection 2 


Lee 


A= 


2 


)xt0=8.5m, P, =7+10 =10.16m 
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Ae YS, 8.5°°V0.001 _ 1.12 j 
= —_V¥0 -§_—) N= = 682m / 5 
nP 0.035x10.16"" 0.164 


Subsection 3 


Q, 


=(4*7).40=5.5m?, PB =4+10=7.16m 
2 3 


_ A/S, 5.5°°J0.001 _ 0.542 


= __V 0 = =~ =417m?/s 
Q; nP? ~—-0.035x 7.1677 0.130 


Hence Q=Q1+Q2+Q3=86.1+6.82+4.17= 97.09 m3/s 


2. a) For the best rectangular section, b=2y,so P=b+2y=2b, A=by= b*/2 


A’? IS, — (b?/2)°2.40.002 


Q=— pe 0.015(2b)22 
So 
b*’? =16.91 b =2.89m 
Channel width is 2.89m and depth is 1.45m 
b) 
A= by =2.89x1.44 = 4.16m* 
If y=1.0, then 


b=A/y=4.16/1=4.16m, P=b+2y=4.16+2=6.16m 


_ WPS). ane? /0,002 


= = =9.55m?/s 
Q np?" 0.015x6.167’° 


3. For the same A, the best hydraulic sections for rectangle, trapezoid and semicircle are 


2 
1) Rectangle A= am so b=V2A, p=2b=2V2A=2.83VA 


2) Trapezoid A=(b+xy)y= = + Say = say" 


V3y v3 


So y= a =0.76VA 
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So semicircle has the least perimeter among all sections with the same area; hence it is the most 
hydraulically efficient of all sections. 


4. Method 1 
The area and wetted perimeter of a triangle are 


A=)", 60 y= Asx 


P=2yy1l+x? =2VAVx +1/x 


oP a1" = 


dx V1/x+x 
ie, 1-1/x* =0 
x=1=tan 45° 


This means that the section is a half square. 
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Method 2 
A=y’ tan(a@/2), y=. /A/tan(@/2) 


P=2y/cos(a@/2) = 2,/A/tan(a@/2)/cos(@/ 2) 


To minimize P is equivalent to minimize P’, i.e., 
Min z =1/[tan(@/2) cos’ (a@/2)]=sin ‘(a@/2)cos '(a@/2) 


dz/da =[-sin *(a/2)cos(a/2)cos ‘(a/2)]/2+sin ‘(a/2)cos*(a/2)sin(a@/2)]/2 =0 
—sin *(a@/2)+cos *(a/2)=0 
tan(a/2)=1 so a@/2=45° hence a =90° 


This means that the section is a half square. 
Method 3 Using Matlab to do the minimisation 


With P=2y/cos(@/2) = 2,/A/tan(a@/2) /cos(a/2) 


> >fminunc(inline(' 2*(1/tan(x/2)0.5/cos(x/2))'),1) 
1.5709 (radian) i.e., 1.5709*180/3.14159=90° 


Or with 


P =2yV1+x? =2JVAVx+1/x 
> >fminunc(inline(' 2*(x+ 1/x)40.5'),1) 
Answer =1, hence it is a half rectangular. 
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10. Critical Flow 


Similar to airplanes with supersonic and subsonic speed classifications, water flow in open channels 
can also be classified into supercritical and subcritical flows. As a counterpart to the sound speed, the 
wave speed (also called celerity) is used in hydraulics. 


Push 


10.1 Small Wave in Open Channel ) 


A small wave can be generated in an open channel by 
pushing a plate. The celerity of a small wave in open 
channel is 


c=+v9Dm 
The derivation of this formula is explained in the unsteady flow chapter. 


If a pebble is dropped into a channel with flowing water, the pattern generated could be one of the 
following cases that represent the four types of flow. 


Still water Critical flow Subcritical flow Supecritical flow 


10.2 Critical Flow 


Critical velocity 


ae 
gD, 
where 
Vv the mean velocity 


Dm __ hydraulic mean depth, A/B 
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Flow state 
Critical Fr=1 or V=JgD,, 
Supercritical Fr>1 or V >JgD,, 
Subcritical Fr<1 orV <gD,, 


In relation to critical flow 


Ye the critical depth 


Y>Y¢, Subcritical Y<Y» Supercritical 
Sc the critical slope 

So<Sc, mild slope So>Sc, steep slope 
Ve the critical velocity 

V<Vc, tranquil flow V>Vce, rapid flow 
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10. 3 Critical Flow Computation 
a) If the critical depth is known, V¢ is solved by 
V.=/gD,, and Q=AVc 


b) If the discharge is known, yc is solved by (trial and error) 


Q=A,gD,, 


c) To get the critical slope Sc , use the Manning Equation 


R 2/3 iS 
V.=—— 
n 
2 
. . nV. 
Rearrange it into S.= R23 
Example 


For a trapezoidal channel with a base width b=6.0 m, side slope x=2 (i.e., 1 vertical : 2 horizontal ) and 
Manning’s n=0.02, calculate the critical depth, critical velocity and critical slope if its discharge 
Q=17m"/s. 


Solution 


A=(b+xy,)y, =(6+2y.)¥, 


B=b+2xy, =6+4y, 


2h (St2yoye OPIN: 


Ji 6+ 4y. 3+ 2y, 


The discharge is 


Q = AY. = A, V gD, 


3+ 
So 17 =(6+4+2y.)y. gery We 
a+ 2y, 
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Simplifying, 


3 
7.36 = B+ yee 


(3+2y,) 
By trial and error, the critical depth is approximately 


y, =0.84m 


and the corresponding critical velocity is 


A=6.45m* 
Vi= we 2.6m/s 
6.45 


The critical slope 


P=b+2y,v1+x° =9.76m 


R as = 0.66m 
P 


2 2 
s,=( 2%] =| 92°4 ) 0.005 
R 0.758 


Questions 10 
Critical Flow 


1. A rectangular channel is 3.0 m wide, has a 0.01 slope, flow rate of 5.3 m?/s, and n=0.011. Find its 


normal depth yy and critical depth yc. 


(Answer: 0.4m, 0.683m) 


2. Water flows in a rectangular channel at a depth of 1.22m and discharge of Q=5.66m°/s. Determine 


the minimum channel width if the channel is to be subcritical. 


3. A rectangular channel has a bottom width of 8 m and Manning’s n=0.025 
a) Determine the normal slope at a normal depth of 2m when the discharge is 12 m’/s; 
b) Determine the critical slope and the critical depth when the discharge is 12m*/s; 


(Answer: 1.34m) 
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c) Determine the critical slope at the critical depth of 1.5m and compute the corresponding discharge. 
(Answer: a) 0.0002, b) 0.009, 0.61 m, c) 0.008, 46.1m*/s) 


4. For a trapezoidal channel with a base width b=3.0 m, Manning’s n=0.025 and side slope x=2 (i.e. 1 
vertical : 2 horizontal ), calculate the critical depth, critical velocity and critical slope if its discharge 
Q=10m’/s. 

(Answer: 0.86 m, 2.46 m/s, 0.0076) 
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Solutions 10 
Critical Flow 
1. From the Manning equation 


A= by=3y, P=b+2y=3+42y 


5/3 
Since o- Ay. Xe) (3y)*" [0.01 
53= : 


nPp23 3= 
0.011(3+2y)*” 
Simplifying 
5/3 
0.093 = _by 
(34+2y) 
Trial and error 
y=1.0 0.342 
y=0.5 0.125 
y=0.3 0.057 
y=0.4 0.089 
y=0.45 0.107 close to 0.4 


So yn=0.4 m 


From the Froude number=1 


Q =AJgDn = By W9c 
Q 
3/2 __ 


Ve ae ox 
BYa 


1 
5.3°\3 
Soy, -( = 0.683m 
3 q) 


2. For a critical flow 


Fr= =‘ 
Joy 
5.66 / (1.22b) -1 


Vl22xg 


b=1.34m 
Hence, the minimum channel width if the flow is to be subcritical is 1.34m. 


3. a) For the given data 
A= by =8x2=16m* 
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P=b+2y=84+2x2=12m 


From the Manning equation 


2 2 
0.02512 
S. | og = ses = 0.0002 
AR3) — \46x1.333 


Thus, So=0.0002 will maintain a uniform flow of 12m’/s in this channel at a depth of 2m. 


b) The critical depth 
A= bye = Bye 
A by, 
D ow = per dies a 
mB BH ve 
Flow velocity is 
Q 12 15 


and at the critical depth 


Vo= gDn 
1.5 

So —_ =VIVc 
Cc 


1.52 
Yo= ye —0.61m 
g 


Then A= by, =8x0.61= 4.88m" 
P=b+2y, =8+2x0.61=9.22m 


2 2 
Qn 120.025 
>| =|——s| = 0.009 
ARG 4.880.533 


c) A=by,=8x1.5=12m? 
P=b+2y, =8+2x1.5=11m 
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and at critical flow 


v= GD = VI = J1.5g = 3.84m/s 


The critical slope will be 


2 
Von 
S. = | = 
R3 


The flow rate is given by 


2 
3.84 x0.025 
—s— | =0.008 


1.095 


Q = V,A=3.84 x12 = 46.1m? /s 


A= (D+ XV e)V¥o = (3+ 26) Ve 
B=b+2xy,=3+A4y, 

A _ Gt 2c 

B 3+ 4y, 


The flow velocity is 


Dae 


10 
oe 


A (3+2y)¥e 


and at the critical depth 


Vo = 9D in 


10 _ [,G#2y Ive 
(3+ 2Yc)Ve 3+ 4Y¢ 
ea 
Simplifying, 10.19 = [e+ 2x del 


(3 +4yc) 


sO 


By trial and error, the critical depth is approximately 


Yc=1.0 17.86 
yc=0.5 1.6 
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Yc=0.75 6.41 

Yc=0.85 9.96 

Yc=0.86 10.38 close to 10.19 
So Yo = 0.86m 


and the corresponding critical velocity is 


A= (3+2 x 0.86) x 0.86 = 


The critical slope 


4.06m? 


P=b+2yWV1+x* =3+2x0.86V1+2° =6.85m 
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11. Rapidly Varied Flow 


For a rapidly varied flow, its water surface profile changes suddenly. This occurs when there is a 
sudden change in the geometry of the channel (to be described in this chapter) or in the regime of the 
flow (hydraulic jump, to be described in the next chapter). 


11.1 Sudden Transitions ——_—_:«. 


Transitions with a change of cross-sectional dimensions occurred 
in a relatively short distance 


The energy equation is used to solve these problems. Se WN \ 


The Energy Equation, 
y2 
H=y+—+z 
2g 


Total energy head = pressure + kinetic+ potential 


Conservation of energy 


(Usually H]gss is ignored in the computation of sudden transitional flow) 
ie. H, =H? 


11.2 Depth of Flow 


A steady uniform flow in a rectangular channel of width 5m is 
interrupted by the presence of a hump of 0.2m in the channel. 
The upstream depth is 1.25 m and the discharge is 10 m’/s, find 
the depth of flow at position 2 (taking a=1 and assume no energy 
loss). 
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Solution: 


Applying the energy equation, 


Known: y;, Vj, Az 
Unknown: y9, V> 


Since 
L= pe 5s 1.6m/s 
by, 5x1.25 
by, 5y, Yo 


So the energy equation would be 


2 2 
ean e: 
2g 2gy 


2 


1.25+ 


Simplify it into 
y; -1.18y; +0.204 =0 


This is a cubic equation, in which y, is the only unknown (with three roots). Solve the equation and 


find yo. 
y= -0.36, 0.59, 0.96 m 


Only one root is the right answer. The negative one can be eliminated first. 
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The critical depth for the channel 


Q= AV, = AJ@D,, = by../ay. 
10 =5y,,/gy. 


so 


Since the critical depth is the limit for the flow over the hump, the root 0.59 m can be eliminated. 


The right answer for y2 is= 0.96m 
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Rapidly Varied Flow 


11.3 Height of Hump 


If the depth of flow in a 2m wide rectangular channel is 2.5 m 
when the discharge is 4 m’/s, estimate the height of the hump 


which would produce critical conditions. 
Solution: 


The energy equation 


2 2 
y; Pe = y. + Veo AZ 
2g 26 
Vas : = 0.8m/s 
A, 2x%2.5 
Q= AV, = AV gy. = bY. VD 
So 
4 1/3 
y, = (=| =0.742m 
g 
2 : =2.70m/s 


© by. 20.742 


Substitute it into the energy equation 


2 2 
2.5+ ne = 0.742 + aly + Az 
2g 2g 
2.53 =1.114 Az 

Then 
Az =1.42m 
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11.4 Specific Energy 


(Pressure + Kinetic energy) 


2 


V 
E, = Vr ya = 
2g Q= constant 


At critical depth, the specific energy is at the minimum point. 


: eee _ 
: ae PEs 
Mathematical derivation: 
EF 2 2 
From - -=0 and E, =yt=yte 
Ly g g dy B 
sO 
1 Ps 2 2 
OE EA ps cg -2 4 (2 )-1-224 
dy dy 2gA 2g dy\A 2gA° dy 
dA 
From the figure on the right, dA = Bdy 
2 2 
: == ia, 1 SR a0 
gA gD, 


Hence, at the critical depth, the specific energy is the minimum. ) 


Questions 11 
Rapidly Varied Flow (Sudden Transitions) 


1. Water is flowing at a normal depth in a 3m wide rectangular channel with a bed slope of 1:500. If 
Manning's n=0.025 and the discharge is 5m°/s. Calculate the height of a hump which would produce 
the critical flow without causing backwater upstream (i.e., raising the upstream water level). 

(Answer: 0.31 m) 
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2. Water is flowing at a velocity of 3.4 m/s and a depth of 3.4 m in a channel of rectangular section 
with a width of 3.4m. Find the changes in depth produced by 


1) A smooth contraction to a width of 3.0m 
2) The greatest allowable contraction for the flow to be possible upstream as described. 
(Answer: 3.0m, 2.89m) 


3. The normal depth of flow in a rectangular channel (2m deep and 5m wide) is 1m. It is laid to a slope 
of 1m/km with a Manning’s n=0.02. Some distance downstream there is a hump of height 0.5m on the 
stream bed. Determine the depth of flow (y:) immediately upstream of the hump and the depth of flow 
(y2) above the hump. If the hump is reduced to 0.1m, what values will y; and y2 be? 

(Answer: 1.25m, 0.54m, 1.0m, 0.87m) 
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Solutions 11 
Rapidly Varied Flow (Sudden Transitions) 


1. From Manning equation 
A= by=3y, P=b+2y=3+2y 


2 5 
A3.f5 3 
since 9 =f2N gg 5 = Gy) *V0.002 
ae 0.025(3 +2y fh 
Simplifying 


2 
2.795 =O), 


(3+2y)3 

Trial and error 

y=1.0 2.13 

y=1.2 2.74 close to 2.795 
So y=1.20 m 

5 
= x = =139m/s 
Ay ax 12 


From Froude number=1 


OQ =A eDm = By Vc 


y?? =Q/(BYg) 


sO 
y. =[5°'/(3°g)I'" =0.66m 
V.= 2 =2.53m/s 
y-B 


From the energy equation 
V7 Ve 
+— =Az+ y+ 
JI 2¢ Ye 2¢ 


2 2 
12+ aaa = Az+0.66+ za 
2x98 2x98 
Az=0.31m 
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2. 1) By the energy equation 


2 2 
y+ — = Yo +e 
2 22 
2 es 
3.4+ =y,+—> 
2 2g 
Since Q = V,y,b =3.4 x 3.4x 3.4 =39.3m’ /s 
39.3 13.1 
y,--2- -393_ Bd 
by, 3Y. Wy 
So 
3.4 13.1° 
3.44 =y, ; 
28 2 8Y> 
Simply 
8.75 
3.99 = y, +> 
a) 


Before trial and error, find the critical depth 


From V.=AJgy. oe 2 2\71/3 2 72 
y.b y, =[Q* /(gb*)]"* =[(39.3* /3° g)]"" = 2.60m 


c 


by trial and error between (2.6 m to 3.4 m) 


y2=3.0 m, 3.97 
y2=3.1 m, 4.01 


So y2=3.0 m 


2) Critical depth at section 2 


2 V 2 
Energy equation ytt=y,t+ > 
2 22 
A v 
3.4+ =y,+ ae 
22 “ 2g 
7 a 
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3.99 = i, y. = 2.66m 


9 
b, = aes = 289m 
V Ne WV 


Maximum contraction to 2.89 m 


3. 

a) Area 

A=5m* 

Hydraulic radius 
A a 


R=—=——=()./Im 
P 542 
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Q=AV=A = 629m? /s 


a ae 0.71" V0.001 
n 0.02 


Q= AY. = (by) gy. =bY gy.” = 6.29 
y, = 3/(6.29/b)* /g =0.54m 
V.= gy, = (054g =2.3m/s 


b) Check the height of hump that can produce critical flow without causing backwater 
The total energy at the hump 


2 2 
yo aoe +z=054+ = +z=081+2z 
2g 2 
Before the hump 
2 2 
y,+ us =1+ (se) = 1.08m 
2g 2g 
So 
081+z=108 


z=108-0.81=0.27m 


c) Therefore, 0.5m hump will produce critical flow over the hump and cause backwater 
The total energy at the hump 


Vv? 2 
y,+—4+2=054+—+05=131m 
2g 2g 
Before the hump 
2 2) 
ip ay Oe. =1.31m 
2g 2g 


y, + 0.081y,* =1.31m 

Trial and error (y;>1.0m due to the hump): 
yi=1.1, LHS=1.17 

yi=1.2, LHS=1.26 

yi=1.3, LHS=1.35 

yi=1.25, LHS=1.30 

So yi=1.25m 


d) If the hump is 0.1m, no backwater will be produced, therefore, 
y, =1m 
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The total energy at the hump 
Vy -2 
Y> vont Z=y,+0.081y,° +0.1 
g 


Before the hump 

2 2 
Vi shee (6.29 /5) 7 
2g 2g 


108m 


Yt 


SO 

es 0.081y,° +0.1=1.08 
ie; 

y3 = 0.98 y; +0.081=0 
Estimate y> (trial above 0.54m) 


y, =0.87m 
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12. Hydraulic J ump 


A hydraulic jump occurs when a supercritical flow meets a subcritical flow. The resulting flow 
transition is rapid, and involves a large energy loss due to turbulence. The depth before the jump is 
called the initial depth and that after the jump is called the sequent depth. 


It is used for dissipating energy in flow to protect dams, weirs, and other hydraulic structures. 
12.1 Jump Equation 


Because significant energy loss in hydraulic jumps, the energy equation is not a suitable tool in 
hydraulic jump analysis, therefore the momentum equation is applied. 


The Momentum Equation 
F = pQVV, -V,) 


In flow direction 


Net force =F1-F2 


(ignore the channel frictional force) 


Momentum change = Mj-M 1 


Hence 
R-R=M)-M 
sO 
F+M,=5E+M, 
Le. 


F + M =constant 
Hydrostatic pressure force 


F=pgAy and M=pQV 
where A = area y= the depth of centroid 


sO 


WAY, — AY.) = QV, — Vi) 


Download free ebooks at bookboon.com 


114 


Concise Hydraulics Hydroaulic J ump 


The jump height can thus be evaluated knowing either upstream or downstream conditions 


Energy dissipation 
2 


V, Vv 
AR =F, =Hy= +2) =O44— 2) 
2g 2g 
Hydraulic jump lengths can only be roughly estimated. One of the popular equations is L=6.1 y2 


Examle * 


A hydraulic jump forms in a rectangular open channel 125m 


with a width of 2 m in which water is flowing at 4 m/s. yi=? 

If the depth just downstream of the jump is 1.25 m. v.. ¥ _ 
Calculate the theoretical depth upstream of the jump. SEN SS 
Estimate the energy loss as a percentage of the initial energy. 


Solution 


From R-KRh=M)-M, 


F, = pgAy, = pg(2y,) ] =9.81py, 


= 1.25 
F, = pgAy, aC am =15.33p 


4x4 _ 8 


Q 
M, = pQV, = p(4)— =p p 
A, 2y, y; 


Q 4x4 
M, = 0OV, = 0(4)— = p—_—_-=6. 
2 = PQV, pan Poa ge p 


Substitute them into the momentum equation 


9.81y, 9 -15.330 =6.4p- 23 
Ji 

Simplify it 
y; —2.22y, +082 =0 


Find y;, (trial and error) 
al LHS 
1 -0.4 
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0.5 -0.165 
0.4 0.004 
So y,= 0.4m 


vo 

HA, = y, += =1.694m 
2g 
vo 

A, = y, + =1.380m 
2g 


Energy loss=1.694 - 1.380 = 0.314 m 


12.2 Hydraulic J ump in Rectangular Channel 


If a jump occurs in a rectangular channel, there is no need to use trail and error approach and the 


unknown depth can be solved by the following equations. 


_ {ye ,2a _ [¥2 24° Yo 
y2 = + y= + 
4 gy, 2 4 gy, 2 
Q 


where q= > 


or in another form: 


Y= (JL+8Ri y= (f1+8F) -1) 


For the last example: Input y, =125 and q=2 
y, = 0.396m 


12.3 Hydraulic J ump in Trapezoidal Channel 
A trapezoidal channel with a 6m bed width and side slope of 1 
1:1.5 discharges water at a rate of 8.5 m’/s. If the normal depth is 


downstream of a hydraulic jump is 1.2 m, what is the flow a 
depth upstream of the jump? 
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Solution: 
Since F-F,=M,-M, 
so PIAY, — AY.) = PQAWV, —V) 
From A=y(6+1.5y) 
aE re 
so A, =y,(6+L5y,) : —SE==aa SF dz 


A, =1.2(6+15%12) = 936m’ I. ém | 
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yy 

Ay = [[6+3(y—z)ledz =| 32? +1.5yz? -2° | =3y? +0.5y* 
0 

So 

AY, =3y° + 0.5y° 

AY) =3x12? +0.5x 12° =5.184m* 


Let V=Q/A 
Thus, 
i 1 
9.81(3y7 + 0.5y, — 5.184) = ast( - —_| 
9.36 y(6+1.5y) 
Simplifying 
7.36 

5.97 —3y? —05y? = > — 

y,(6+15y,) 


By trial and error 
y,=0.2 m 
12.4 Hydraulic jump in Sloping Channel 


In flow direction: Net force =F, sinO+F1-F2 
Gravity force 


f, = = PgL(h, +h) F, == pghi cos@ F, == pghi cos 6 


F, Al 
Soh, = A fi+ew; ~1) : 


then vertical depth y,;=h;/cos®@ and y2=h,/cos® 
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Questions 12 
Hydraulic Jump 


1. Water is flowing at a rate of 10m’/s through a rectangular channel 4 m wide, at a depth of 0.5m. A 
weir downstream causes the water to backup the channel and a hydraulic jump occurs. Find the 


sequent depth and the loss of energy at the jump. 
(Answer: 1.37m, 0.234 m) 


2. Water flows in a rectangular channel at a depth of 30cm and with a velocity of 16m/s. Ifa 
downstream sill forces a hydraulic jump, what will be the depth and velocity downstream of the jump? 
What head loss is produced by the jump? 

(Answer: 3.81m, 1.26m/s, 9.46m) 


3. Starting from the first principles, show that the following equation holds true for a hydraulic jump in 
a trapezoidal channel: 


b 2 3 2 
( aes d )+ pQ = constant 
2 3 (b+xy)y 
where y is the depth of flow, b is the bottom width, Q is the discharge and x is the side slope (1 vertical 
to x horizontal). Draw the force momentum diagram for the following conditions and determine the 
initial depth if the sequent depth is 0.2m: Q=50 l/s, b=0.46m, x=1. 


(Answer: the initial depth is 0.04m) 
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Oo 100 200 300 400 500 600 
F+M 


Solutions 12 


Hydraulic Jump 


1. From the rectangular channel equation 


_ |v 2 
V2 = + 
4 gy, 2 
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Input Y, = 0.5m and q=2.5 


Jz = —=137m 
4 9, 2 4 05g 
: ; 25° 
ee ey 4 5 =05+ 2 > = 1774m 
2g 2g; 2x981x05 
: ° 2.57 
sa ae 2 7 = 137+ 7 = 1540m 
2g 293 2 x 981x137 
Energy loss AH=H 1-H 2=1.774 - 1.540 = 0.234 m 
1 
i= oe a3) 
2. Since gi 0.3g 


Use the rectangular hydraulic jump equation 


yy = 2 (fir err -1)= 93 (fivaxo33' -1}=3.81m 


since V,¥,D=V,yb 


és. Yo a sears 
y>, 3.81 
: : 16° 1.267 
Fy Gs oases : ) =9.46m 
2g 2g 2g 2g 


3. 

F + M =constant 

pg( Ay) + pQV = constant 
A= y(b+ yx) 

So 


Ay= [( + 2x(y—z)]zdz = [ bz” {2+ xyz" - 2xz°/3 | = 
by? /2+xy?/3 
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_ PQ’ 


pg( Ay) + pQV = pg(by* /2+ xy*/3)+ = constant 
(b+xy)y 
y 0.25 
Force + momentum curve 02 
0.05 
F+M =1000[g(0.46y* /2+ y* /3) +——_] 0.15 
(0.46+ y)y 
0.05° = 
= 1000[ g(0.23y’ +0.333y’) + ——_- — 
(0.46+ y)y 0.05 
0 T T 
y=U2 FeM=135;y-0.15 FM =o9; 0 100 200 300 400 500 600 
y=01 F+M=70.5; y=0.05 F+M=104; van 


y=0.01 F+M=532; y=0.03 F+M=172 


So the initial depth is about 0.04m 
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13. Hydraulic Structures 


Hydraulic structures are devices which are used to regulate or measure flow (weirs, flumes, spillways, 
gates, etc). The energy equation is applied to estimate the discharge and the momentum equation is 
used to calculate the forces on the structures. 


13. 1 Flumes (Venturi) 


Flumes are used to measure discharge in channels. The energy equation is used to derive the flume 
formula. 


throat 


exit | 
il 2 
If the width b is small enough, there will be a critical flow in the throat (also it will cause backwater 


upstream to get extra energy and a hydraulic jump after the flume to get rid of the energy gained from 
the backwater). 


The energy equation between section 1 and 2 (ignore energy losses): 


A,=H 


1 pais 


ve (Jav.) 


The total energy headat 2 is H,=y,+——=y,+ =1.5y. 
2g 2g 
H 
So -—1 
a 153 


The discharge through the flume will be 


5 
Qhaeat = AV, = by, VV IY. = bg y2” = ba( =| =1.705bH;° 


In reality (to consider the energy losses) 


0=C.0..,=17050C.4,° 
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where Cd is the coefficient of discharge for the energy losses (see BS3680 4C) 
If the depth h is used Q=1.705bC,C,h'” 


where Cv is the coefficient of velocity for the upstream velocity head (also see BS3680 4C). 

Such a structure is an example of critical depth meter. 

Example Venturi Flume 

An open channel is 2 m wide and of rectangular cross section. A venturi flume having a throat width 
of 1.0m is installed at one point. Estimate the discharge: 

a) if the upstream depth is 1.2 m and a critical flow occurs in the flume 

b) if the upstream depth is 1.2 m and the depth in the throat is 1.05 m 

(Take Cv=1 and Cd=0.95) 


Solution: 


a) If the flow in the throat is critical 
Q= 1.705C,C,b,h'° =1.705x 0.95x1.0x1.2'° =2.13m?/s 


b) If h=1.2 and y,=1.05, use the energy equation (ignore the energy loss) 


y a 
y —=y ee 
one ag 

2 2 

(Gia 295, 

2g(2x1.2) 2g(1x1.05) 
Q? 
(0.907 —0.174) = 0.15 
2g 


So.Q..,=2.0m'/'s 
Considering the coefficient of discharge with energy losses 


Q =0.95x 2.0 =1.90m?/s 
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13.2 Weirs (Broad-crested weir) 


In a similar way to Venturi flume 

Q=1,705C,,C bh” 

where: 
Cd — the coefficient of discharge 
Cv — the coefficient of velocity 


Broad crested weirs are robust structures and usually used in larger rivers in comparison with flumes. 
However, flumes usually create relatively smaller backwater (or afflux) than weirs and hence less 


energy loss. The floating debris can usually pass through them smoothly (a useful feature in sewerage 
treatment works). 


13.3 Energy dissipators 


Kinetic energy in flow can erode the foundation of dams and other hydraulic structures. It is important 
that excessive energy should be dissipated to protect the foundation. 


1) Spillways 


are safety structures for dams. They must have the : 
. . ‘ : Spillway 
capacity to discharge major floods without damage to 


Hydraulic jump 
the dams. 


Q=Cbh’” (C usually lies within 1.6 - 2.3) “ae ee 
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2) Others a) Stilling basin; b) Submerged bucket; c) Ski jump 


Co Be 


Stilling basin calculation 


The velocity at Section 1 


v= p./2gh, 


¢ is usually taken as 0.90-0.95 


The depth before the jump (rectangular channel) 


mS 
bv, 


Mi 
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The sequent depth 


y= + 8F; ~1) 


The wall at the end can be treated as a weir, so from Q =1.705bh"” (ignore the energy loss) to derive 
h. 


The height of the wall will be Az = y, —h 
The length of the basin is shorter than free jumps (6.1y>2) due to the baffling effect from the end 
wall, so 


1, = (4.0 ~ 5.0)y, 


13.4 Sluice Gates 


Energy Equation (Section 1 - 2) 


1 2g 2 2g 


If y, and y, are known, 


discharge through the gate can be derived 


Q=by,y, 


1 V2 


In practice, it is easier to link the discharge with the gate opening, so 


Q=C,by,/2gy, 


where 


C 


Cc 


a a 
Ye elg and tI Cel Y 


Ce can be found from hydraulics manuals for various types of gates (usually between 0.61-0.66) 


The force acting on the gate can be derived from the momentum equation between section 1 and 2 
(ignore frictional forces). If the downstream depth is greater than the critical depth y,, a hydraulic jump 
will occur in the channel. The sequent depth ys; can be derived from the depth at y2. If the 
downstream depth is greater than y3, a submerged gate flow will occur and the flow equation derived 
here is not applicable anymore. 
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Questions 13 
Hydraulic Structures 
1. The depth-measurement scale of a Venturi flume with a 1.2m throat give an h reading that is 0.02m 


too large. Compute the percentage errors in discharge when the observed h readings are 0.2m and 
0.6m (assume the critical flow condition for both cases). 


(Answer: Around 17%, 5%) 


2. Determine the width bp of the throat for a Venturi flume so that it will just produce critical flow, 
given Q=40m°/s, bj=10m and h,=2m. Also calculate the possible downstream supercritical flow depth 
hs (Neglect energy losses). 


(Answer: 7.17m, 0.75m) 


6007 'A'G Swewsds Wy] 282] @ 


iB \ 
3 It’s only an 


opportunity if 
you act on it 


IKEA.SE/STUDENT 


. 


ven oes a : 
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Solutions 13 
Hydraulic Structures 


1. Method a) 


From 
Q=1.705C,C.b,h°” 
and 
d 3 
~ = 1.705C,C,by 5h"? = 2,558C,C,b,h"” 


then 
312 _ 1/2 
AQe 1.705C,,C.b, 5 ho" Ah= 2.998C,C,b,h Ah 


The percentage error will be 


AQ _ 2.558C,C,b,h"?Ah _ 2.558x0.02 _ 0.03 


Q,, 1.705C,C.bh®? -*1.705h_ ~~ 
At h=0.2 m 

ae = ON =16.6% 

Qos 0.18 

At h=0.6 m 

AQ 0.03 5 nog 

Qc 9.6—0.02 

Method b) 

From 


Q=1.705C,C,b,h*” 


If h is 0.02 too large, the percentage error will be 


AQ _ 1.705C,C,b,(h)'* -1.705C,C,b,(h—0.02)'5 _ (hy'® -(h—-0.02)!° 


A Ine aA A moanst5 moanstS 
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AQ _ 0.2'° —0.18'° 


“—__ =17.1% 
Q 0.18" 
At h=0.6 m 
15 1.5 
AQ _ 0.6 ee _5.9% 
oe 0.58 , 


2. From the Venturi flume equation 


Q =1.705b,H;"’ 


vy Oo 40 

Since Q=40 m3/s and H, =h,+—-=2+ 5 =2t say = 2.204m 
2g 2g(b,h) 2g x10° x2 

So Q = = = 7.170m 


b, = 
*-1.705H?? = 1.705 2.205'° 


(You can also derive by from the energy equation) 


For the downstream depth h3, 


‘ 8157 

fais ek eo 
2g 29(b,h,) h; 
0.8157 


Since H,=H, So 2.204=h, + 2 


3 
Trial and error (below 2m), 


h, =0.75m 
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14. Gradually Varied Flow 


Gradually varied flow refers to the flow whose depth and velocity smoothly change over a long 
distance. 


14.1 Equation of Gradually Varied Flow 


The energy line slope S¢ 
dH 

Ss; =-—— 
dx 


( The ‘-’ sign is to compensate for the negative 
dH so S¢ is positive) 


The channel bed slope is Sg 


s,-—# 
dx 
Th aoe ( — +2) 
e energy equation ce = de y 2g 


Only the kinetic energy change needs special 
treatment (other terms are straightforward). 
d 
ee ae 
dx dx 2g 


For the kinetic energy change 
d ae d (Q/ Ay", _Q dA~ 
dx 2g dx 2g 2g dx 
Q’ dA Q° Bdy 
~ gA® dx _ gA° dx 


(Since dA=dy B) 


5 _Y_ QBdy_ 
= ' dx gA’dx ° 
ve *B 
Since Fr* = ms . 
gD,, gA 


hence, the final general flow profile equation is 


dy _5,-5, 
dx 1-Fr’ 
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From the Manning Equation (Replace S, with S; and S¢ # So for gradually varied flow) 


2/3 | 
ve R S; a nv? 
~ n sO f R?? 


The qualitative relationship between the relevant variables is as below. 


ml tl a 
dx 1-Fr’ 
2/3 
n gD,, 
When y> yn Sf < So When y> Yc Fr<1 
When y< yn Sf >So When yK< Ye Fr>1 


. iy 
Pa 


4 


fe} 
RURAN CE 
ALAIN 
NORLD 
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14.2 Classification of Surface Profiles 
Flow surface profiles are classified based on the actual depth and the channel bed slope. 
1) Depth 


The depth space is divided into three regions by normal and critical depths. 


region 1 
i ee region 2 Yn 
region ) ets Ye 


ee 


2) Slopes 


The channel bed slops are classified into five groups: M mild, S steep, C critical, H horizontal and A 


adverse. 
Region 1 Region 1 Region 1 
i eee ie ee Ye 
SeSSecn Region 2 Yn WRT ea Region 2 Ye 
Re Region 3 Yn 
Region3 7 Ye Region 3 Yn 


Critical slope 


Mild slope Steep slope 
Region 2 Region 2 
Ueee-or Yc 
sensecsscecasesecesscnesecenesees Ye a eeee 
: Region 3 
Region 3 


Horizontal Slope Adverse slope 


(No normal depth) (No normal depth) 


3) Profile name 


A letter+ adigit eg. M1, S2 
where: Letter - the type of slope (M, S, C, H or A Digit - the region (1, 2 or 3) 


14.3 Flow Profile Sketch 


It is possible to sketch flow profiles based on the analysis of the general flow profile equation. Here is 


an example for a mild slop case. 
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Example: In a mild slope, if y is in Region 1, which profile is correct? 


dy 


= , 
i ie Yo M 
It is subcritical flow, so Fr<1 


-- 
~eae 
-. 
_ 
~-LL 


and y>yp, so Sf < Sg 


S -—S ee 
dy _So~S; _ positive | 9 ee 


dx 1-Fr’ positive 


The profile with increasing depth is the correct answer. 


Self Practice 
For the same channel bed slope, if y is in Region 2, which profile curve is correct? 


-~-. 
7 see 
-. -- 
ee - 
- -- 
asia --2 


(Answer: the one with decreasing depth) 
Steep slope cases (derived in the similar way as mild slope cases) 
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14.4 Flow transitions 


1) A smooth surface link 
Subcritical to Subcritical 
subcritical Subcritical to Supercritical 


Supercritical to Supercritical 


2) Linked by hydraulic jump 
Supercritical to Subcritical 


supercritical 


subcritical 


S.<S¢ 


Questions 14 
Gradually Varied Flow 


1. Water flows in a rectangular channel with a bottom slope of 0.8m/km and a head loss of 0.4m/km. 
At a section where the flow depth is 1.7m and the average velocity 1.8m/s, does the water depth 
increase or decrease in the direction of flow? 
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2. Roughly sketch the possible flow profiles in the channels shown below 


a) 


WHAT'S MISSING IN THIS EQUATION? 


t 

S = 

xe) ~ : ‘or 

S i “4, Ps > 

oO : <—/ ~ P > - 

= AS - “< < 

c | as 3 fate ’ 
é a 

3 See FD 

Oo Dre <2: = 

; : 

wn 

oO 

a 

a 
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b) 


Steep Mild 
Mild 


3. Show that the gradually varied flow equation is reduced to a uniform flow formula if dy/dx=0 


Solutions 14 
Gradually varied flow 
1. The Froude number 


re Vv _ 18 
Jay Jl.7g 


From the general flow profile equation 


= 0.44, 5,=0.0008, S;=0.0004 


ho ie = 

dy ac) Pa f 7 0.0008 ee ~0.0005>0 
dx l= Fr 1—0.44 
so the flow depth is increasing. 
2.a) 

a ae T 

Mild 
Mild 


Steep 
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b) 
Steep Mild 
Mild 
dy 5,75; 
3. From dx. 1-Fr 
cn 
If Oy dies (a S,-S,=0 
x =F ; 
nQP??\" 
Since S,= ~e 
nQP2!?\" 
So a -( AD ) 


Hence Tt ic the Mannino Fanation far nnifarm flawe 
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15. Computation of Flow Profile 


15.1 Introduction 


dy _ 5, =S¢ 
dx 1-Fr’ 


The general flow profile equation is 


It is possible to derive y=f(x) from a simple differential equation. For example, if dy/dx=0.5, the 
solution for it would be y= x/2 +c, where c is a constant. There are an infinite number of curves by 


changing the c value. To uniquely identify the profile, a control point is required (i.e., a boundary 


Cy y=x/2+ 3 

ee Co ee 

a (0,3) 
C3 

a PSs ated Point 


SSS xX 


condition). 


For example, if y=3 when x=0, the curve function would be fixed as y= x/2 +3. 
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For practical flow profile computation, there are three types of solutions: 1. direct integration; 2. 
graphical integration; 3. numerical integration. The first two methods are superseded by numerical 
integration methods due to the widespread use of computers. 


15.2 Numerical integration methods 


1. The direct step method: _ to solve section length from depth. 
2. The standard step method: to solve depth from section length. 


1) Thedirect step method = 
Known: y, and y, Unknown: AX=? fe Jay 


By using a finite difference form 


1—Fr? - 
ax-ay[2 ae | 
0 f mean 


where 'mean' refers to the mean value for the interval (Ax). 


2) The standard step method 
Known: y, and Ax Unknown: y,=? 
By using a finite difference form 


-S 


Ay = ax{ 2 f (Iterations are required since y» is unknown and the means of S; and Fr need 
mean 


1- 


y2 +) 
For approximation, use 


ae 
Ay = Ax 1— Fr to avoid iterations. 
Section1 


15.3 Computation Procedure through an Example 


Determine the flow profile 
for the following conditions 
if 10 depth steps are used in 
the calculation: 


Channel: Q=600m”/s, 
rectangular, width 50 m, 
So=2m/km, n=0.04 
Weir: CdCv=0.88 

Sill height: 2.5 m 
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Solution: 


a) The normal depth yy from Manning's equation 


ANS, 


e. prem y'°V0.002 


sO 273 


nP* 0.04 x (50 + 2y, ) 


By trial anderror y, =4.44m 


b) The critical depth yc 
2 2 1/3 
Q = AV, =by../gy. » sO ya po ={ 20/80 = 2.448m 
g g 


c) The depth over the weir (h) 
From the weir equation Q=1.705C,C,Bh*” then 600=1.705x0.88x50h"* so h=4m 


d) The total depth over the weir (hy) 


Yw=h+2.5= 4+2.5=6.5 m 
Yw>yn>yc. This is a M1 type profile. 


e) Profile computation 


Take yw=6.5 m as the initial depth and y=4.5 m (slightly greater than yy) as the final depth, and 


proceed upstream at a small interval of depth 
Ay = (4.5-6.5)/10 = - 0.2 m 


Simplify the formulas to be used as 


V Q 600 3.831 


Joy byslay ~ S0yVay - yy 


Fr = 
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nQP?"\ (0.04x600P°* \ 576p*™™ 
S; ~ Ae? ~ At = Al 


1- Fr’ 
Ax = Ay 
S,-S 


zs] 


0.2312 


0.2423 


nN 


Dud 0.2543 


9.8 


59.4 


Note: the negative values of x indicate that they are 


in the opposite direction of X. 


15.4 Further Computational 
Information 


-2000 -1000 0 


1) Computation steps for flow profiles 

a) find the critical and normal depths; 

b) find the slope type; 

c) find a control point with a known relationship between depth and discharge (weirs, flumes 
and gates): 

i) Subcritical flow: its profile is controlled from downstream; 


i) Sunercritical flaw: its nrofile is cantralled fram instream 
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2) Computation direction 


Computation direction 


Yc 


Mild slope 


Upstream if the flow is subcritical 


Control point 


Downstream if the flow is supercritical 


Computation direction 
Ye 


Steep slope 
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3) Locations of the control points 


Subcritical flow 


Control point Control point 


Mild slope y = 


yw is calculated from the weir equation (plus weir sill 


Supercritical flow 


Control point 


ge eee 


Steep slope 


Steep slope 


a 
~~, Yc 


Questions 15 
Computation of Flow Profile 


1. A rectangular channel is 3.0 m wide, has a 0.01 slope, discharge of 5.3 m°/s, and n=0.011. Find yn 
and yc. If the actual depth of flow is 1.7 m, what type of profile exists? 
(Answer: 0.4 m, 0.683 m, S1 curve) 


2. A rectangular channel has a gradient of 2 in 1000 and ends in a free outfall. At a discharge of 
4.53m°/s with Manning's n=0.012 and b=1.83m, how far from the outfall is the depth equal to 99% of 
normal depth if five depth steps are used in the calculation? 
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Solutions 15 
Computation of Flow Profile 


1. From the Manning equation 


A= by=3y, P=b+2y=342y 


3 
5 
Aaa 3y)z V0.01 

Since Q = VS so ee 

ss 0.0113 + 2y)3 
Simplifying 

5/3 
0.093 - 


93 = —+"—_ 
G42y)" 


Trial and error 


y=1.0 0.342 
y=0.5 0.125 
y=0.3 0.057 
y=0.4 0.089 
y=0.45 0.107 close to 0.4 


So yn=0.4 m 


From Froude number=1 


Q=AYgDm = BycW Ve 


Since y=1.7 m,S0 Vy <Yo<y 
This is a steep slope and the profile is S1. 


2. From the Manning equation 
5 
A3JS) 


2 
nP3 


Q= 


Since A=by=1.83y, P=b+2y=1.83+2y 


The normal depth 


Download free ebooks at bookboon.com 


145 


Please click the advert 


Concise Hydraulics 


Computation of Flow Profile 


1 
4.53 = —.83y4[ 


2/3 
wae 0.002" * 


0.012 1.83 + 2y, 
Simplifying 
5/3 
0.444 = —_*2___- 
(1.83 + 2y, 
The right hand side 
y=1.0 0.408 
y=1.1 0.463 
y=1.05 0.435 
y=1.07 0.446 close to 0.444 
so yn=1.07 m 


The critical depth 


4.53 = 
Vo = [485 =0.855m 
1.83°g 
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Since yn>yc, the flow is subcritical and the water surface profile is M2. 
The final depth y=0.99y, = 0.99 x1.07 =1.06m 


Ay=(1.06-0.855)/5=0.041m 


Q 4.53 0.79 


i by Jay 1.83yJgy y" 


1.667 3.333 


2 2 
‘ [uae [oo “45g 0.00296 p! 333 
2 ee | 


Al 3 


A A 


0.855 | 1.565 “0.999 0.003588 0 


0.937 1715 370 0871 — 0.002805 9.82 
0.2869 0.000658 
0.978 1.790 3. 0817 0.002511 27.70 


037130 | 


1.019 $1865 33. 0.768 0.002252 - 67.6 


- 0.4430 


106 1.940 (395 (0.724 | 0.002029 /-196.4 


=-196.4m 
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16. Unsteady Flow 


16.1 Basic Types 


If the flow velocity at a point varies with time, the flow is called unsteady flow. They can be further 


divided into (in a similar way as to steady flow types): 


a) Rapidly varied unsteady flows: Surges and bores by dam break, quick operation of control 
structures, tidal effects 


b) Gradually varied unsteady flows: Translatory waves by flood waves, slow operation of control 


structures (gates, etc) 
16.2 Rapidly Varied Unsteady Flow 


1) Dynamic Wave Propagation 


Generation of a solitary wave in a rectangular An observer following the wave 
channel 


The energy equation 


Vv? 2 
+= cea 
MY 2g Yo 2g 


Since yi=y, yo=yth, Vi=c, V2 can be derived from the continuity equation Viy =V2(y+h) 


c cf_y \ 
Therefore y+—=yth+—| —— 
2g 2g\yt h) 
2 hy 
Solving for c, c= gy +h) 
2y+h 
for waves of small height c=q ay 
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Similarly, for non-rectangular channels c=./ gD, 


where Dy - the hydraulic mean depth (A/B) 


2) Surges & Bores 


Vsrg jams 


V2 


Vi —— 
y2 


A surge wave results from a sudden change in flow (e.g., a partial or complete closure of a gate that 
increases the depth), called positive surge (similar to water hammers). It is termed as a bore if caused 
by tides. The momentum equation is used in a similar way to hydraulic jumps. A surge can be 
considered as a steady flow case if an observer is travelling at the surge speed. 
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For a rectangular channel, the continuity equation is 
(V, HV cig V4 = (V, Vs )V> 


or simplified as V_ = me 
v2 7M 


srg 


The momentum equation 
F.+M,=F,+M, (where M = pQV and F = pgAy ) 


There are two equations so the two unknowns can be solved (with the three knowns). Iterations are 
usually required in the nonlinear calculations (assume V.,,.=0, calculate the other unknown from the 
continuity equation, then get Vsrg from the momentum equation, enter it into the continuity equation 
again to get the other unknown, repeat the process until the required accuracy is achieved). Note, the 
surge velocity is much faster than small solitary waves. 


16.3 Gradually Varied Unsteady Flows (Saint-Venant equations) 
1) The Continuity Equation 


At any time t, the outgoing velocity is 


0Q 


— dx 
OX 


Q, =Q, + 
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The net flow into the control volume ( dx long) during dt period is 


(Q,-Q,)dt = [a Jat 
Ox 
The storage change rate (the top part) 
B oy dxdt 
ot 


Applying the conservation of mass 


20, p2Y 4 gy 22, 9A_ 
Ox ot Ox ot 


The equation indicates that when the river level is going up with time (Oy / Ot > 0), the incoming 
discharge is greater than the outgoing discharge (OQ / 0x <0). This represents the flood wave’s 


arrival. 


2) The Momentum Equation (Dynamic equation) 

From the momentum principle, 

External Forces = Momentum Change Rate 

The forces includes the gravity F,, friction F;along the channel bed and wall, and the pressure force Fp 


Gravity 
The total weight of water in the control volume is 
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So 
F, =Wsin @ = pgAdxsin 6 = pgAdxS, 


(when the slope is mild, S, = tan @ = sin @) 


Friction 
The work done by friction force along distance dx is 


F. dx 


This will result in the energy loss of water in the control volume (hj) 
Wh, = pgAh, 
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So 
F.dx=Wh, i.e. F, =Wh, /dx=WS, = pgAdxS , 


S; is the friction slope (use the Manning, Chezy, etc to estimate). 


Pressure 
The hydrostatic pressure force difference between Sections 1 and 2 is mainly caused by the extra depth 
increase at Section 2 to exert extra uniform pressure over cross section area A 


0 
F,= (pgdy) A= [pg Pac] 


The momentum change rate 


O(PQV) 4. 
Ox 


Output - input 


in the control volume with time 


AL(pAdx)V] _ : O(AV) dx 
Ot Ot 


é A(QV)  _aAV 
So pgAS, — pgAS; — pgA 7 =p ( dap oe 
Ax ax at 


Simplify and ignore wind effect 
ayy _ (QV) , AAV) 
Ox Ox ot 


gA(S, -S; 


The right hand side can be expanded as 


A(QV) , AAV) _ [AV 1,62, ,OV_,,AA_ dV, eV 


Q—+V—+A—+V— =Q— +A 
Ox ot Ox Ox ot ot Ox ot 
f ae . 2Q + OA =0 
(from continuity equation Dx ot ) 
So 
dy. eV av 

5,-5 =V + 
Dee a oe 
Rearranging 
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1 ov 
g ot 


_s oy _VOV 
g Ox 


or 


Steady, uniform 


Gravity Friction : With position With time 


Steady, nonuniform Pressure 


Forces Momentum change rate 


Unsteady, nonuniform ; 


The continuity and dynamic equations are also called the Saint-Venant equations 
Solutions 


Unknowns: V, y. Known: Initial and boundary conditions 
Three methods: 

a) Finite differences: currently the most common method. 

b) Characteristics: suitable for the simple regular channel 

c) Finite elements: suitable for irregular boundaries 


16.4 Software packages 
1) Software packages for 1D river flows 


Free Software 

FLDWAV is a generalized flood routing model that can be used by hydrologists/ engineers for real-time 
flood forecasting of dam-break floods and/or natural floods, dam-breach flood analysis for sunny-day 
piping or overtopping associated with the PMF flood, floodplain inundation mapping for contingency 
dam-break flood planning, and design of waterway improvements. Its predecessors are DAMBRK, 
BREACH, DWOPER/NETWORK models that are also available from the same web site. Availability: 
FLDWAYV version 2, June 2000. 

http://www.nws.noaa.gov/oh/hrl/rvrmech/fldwav_instructions2.htm 


HEC-RAS is (formerly HEC-2) is a comprehensive computer program designed to compute water 
surface profiles under a variety of conditions. The HEC-RAS software was developed at the 
Hydrologic Engineering Center (HEC). The system is capable of performing Steady Flow water 
surface profile calculations, and unsteady Flow, Sediment Transport, and several hydraulic design 
computations. HEC-RAS is the most widely used program by practitioners for the calculation of open 
channel flows. Its success is due to the detail that the program is able to handle. Availability: HEC- 
RAS version 4.0, November 2006, http://www.hec.usace.army.mil/software/hec-ras/ 
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Commercial Software 

MIKE 11 is a professional engineering software package for the simulation of flows, water quality and 
sediment transport in estuaries, rivers, irrigation systems, channels and other water bodies. It is a 
dynamic, one-dimensional modelling tool for the detailed design, management and operation of both 
simple and complex river and channel systems. Availability: Purchase from DHI (Danish Hydraulic 
Institute), http://www.dhisoftware.com/mike11/ 


ISIS models open channel systems containing loops, branches, floodplain conveyance and storage 
incorporates standard equations and modelling techniques for structures including weirs, sluices, 
bridges, culverts, pumps, siphons, orifices and outfalls, logical control of moving structures. 
Availability: Purchase from Wallingford Software Ltd, 

http://www. wallingfordsoftware.com/products/isis.asp 


2) Software packages for 2D river flows 


Free packages 

SSITM: SSIIM is an abbreviation for Sediment Simulation In Intakes with Multiblock option. The 
program is designed to be used in teaching and research for hydraulic/river/sedimentation engineering. 
It solves the Navier-Stokes equations using the control volume method with the SIMPLE algorithm 
and the k-epsilon turbulence model. http://www.bygg.ntnu.no/~nilsol/ssiimwin/ 
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Commercial packages 

MIKE21: MIKE 21 is the preferred 2D engineering modelling tool for rivers, estuaries and coastal 
waters. MIKE 21 consists of more than twenty modules covering the following areas: Coastal 
hydrodynamics; Environmental hydraulics; Sediment Processes; Wave processes, 

http://www. dhisoftware.com/mike21/ 


3) Software packages for 3D river flows 


Free packages 
SSITM: as above 


Commercial packages 

MIKE 3: MIKE 3 is applicable for simulations of hydrodynamics, water quality and sediment transport 
in all water bodies where 3D effects are important. MIKE 3 is compatible with MIKE 21 and other 
DHI Software products, allowing for easy model setup, transfer of boundary data, etc. 
http://www.dhisoftware.com/mike3/ 


DELFT-3D : DELFTSD is a fully-integrated two or three-dimensional compound modelling system. It 
simulates the flows, waves, sediments, morphological developments and water quality aspects. The 
sediment transport module models bottom and suspended transport of sediment separately using a variety 
of formulae. The effects of wave motion on transport magnitude and direction are included. The 
morphological module computes bottom changes due to transport gradients and various types of 
boundary conditions. It can be run in a time-dependent way (coupling of hydrodynamics with 
computed bottom changes) or in a time-independent mode. 

In the time dependent mode, animation’s can be made of the bottom development over several years. 
http://www.alkyon.nl/Tools/Delft3D.htm 


Questions 16 


Unsteady Flow 


1. The average depth of oceans is 3720m. Estimate the tsunami wave speed in km/hour based on the 
dynamic wave propagation equation and the average ocean depth. 
(Answer: 687 km/hour) 


2. There is a uniform flow of 4 m°/s of 0.6m deep in a rectangular channel of 4m wide. If the 
upstream gate is suddenly lifted up and the flow in the channel is increased to 6.7m’/s. Estimate the 
surge wave speed and wave height. 

(Answer: 4.5m/s, 0.15m) 
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Solutions 16 
Unsteady Flow 


1. From the dynamic wave propagation equation 


c =,/gy = V9.81x3720 =191m/s = 687km/ hour 


2. The surge wave moves from the upstream to downstream, 


Vorg <=> 
Vi 


V2 


V, = Q, /(by,) = 6.7 /(4y,) =1.68/ y, 
V, = Q, (by) =4/(4x 0.6) =1.67m/s 


The surge can be considered as a steady flow case if an observer is travelling at the surge speed. 


The continuity equation is 
(V, = WA = (V, = Vig > 
_ Vy, -Vy, _1.68-1.67x0.6 0.68 


or simplified as V., = = 


"T (y,-Y>) y, -0.6 y, -0.6 
The momentum equation 


F.+M,=F,+M, (where M = pQV and F = pgAy ) 


-, 1 1 
1:05; 3 Paby: + pQ, (V, =Van) = 5 Paby2 PQ, (V, Vis 


so 2gy; +6.7(1.68/ y, -V,,,) = 2g(0.6)* + 4(1.67 -V,,,) 
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Simplify 19.6y; +11.3/ y, —2.7V,,, =13.7 


srg 
1.84 
Substitute V,, 19.6y; +11.3/y, -—— =13.7 


1 


With trial and error 

y,= 1, LHS= 26, 

0.7, LHS=7.3 

0.75, LHS=13.8 close to 13.7 

so y, =0.75m, the wave height is 0.75-0.6= 0.15m 


What do you see? 


A skateboarder? 


We see an All Options Junior Trader leaving 
the office after a successful day. 
He has the freedom to make his own decisions. 
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17. Hydraulic Machinery 


Hydraulic machinery refers to a device either for converting the energy held by a fluid into mechanical 
energy (turbines) or vice versa (pumps). 


17.1 Hydropower and Pumping Station 


1) Energy of water in motion 


The gross energy head = potential energy head + pressure energy+ kinetic energy 


2 
H=z+ ae + ~ (unit: metre) 
Pg <9g 


Total water power from the gross energy head if expressed in Watt 
P=pgHQ (Wait) 

where 

P= total water power by gross energy head (Watt) 

H= gross energy head (m) 

Q= discharge (m’/s) 

p=density of water (1000 kg/m*) 


2) Efficiency 


Efficiency of a hydraulic machinery device is measured as the ratio of energy output to input. The 
overall efficiency of a hydro power plant/pumping station is a product of the efficiencies of its several 
elements. 


Ns—NyNo...Mn 


Hydropower plant/pumping station output power 


Po = nsPi 
where 
P; = Input power (Watt or K Watt) 


P= Output power (Watt or K Watt) 
Ys= Overall Efficiency 
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3) Hydroelectric power 


Water is dammed and then diverted through a mechanical device to convert water’s kinetic energy into 
rotational energy which can then be converted into electrical energy in a generator. 


Advantages of Hydro Power Limitations 

a) Continuous low-cost production; a) High initial cost; 

b) No consumption of fossil fuel; b) Loss of land suitable for agriculture; 

c) Low maintenance cost; c) Relocation of highways, railroads, even small 
d) No air pollution; towns; 


e) Reservoir can be used for recreation and other | d) Change of local environment; 
purposes (irrigation, water supply, etc.); e) Long transmission line. 
f) Low insurance and low tax. 


System Components 
a) Hydraulic Works; b) Power House; c) Turbines; d) Generators; d) Power lines. 


Classification of Hydropower Plant 

a) Storage Plant 

has a reservoir of sufficient size to develop a firm => 
flow substantially more than the minimum 

natural flow. 


b) Run-of-River Plant 

can use water only as it comes. It is cheaper than 

the storage plant of equal capacity, but suffers as 
seasonal variation of output. 


/ 


c) Pumped Storage Plant 
generates energy at peak load and at off peak, 
water is pumped from the lower pool to upper 


pool. N 
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4) Pumping Station 


Pumps play important role in various civil engineering projects. In water supply, pumps are necessary 
if gravitational flow could be not achieved. They are also used in urban sewers, drainage of low land, 


abstraction of water from bore hole, etc. 


Upstream storage reservoir 


Delivery pipe 


GUE 


Suction pipe 


Screen 


A pumping station 


Classification of pumping stations 

a) Abstraction from surface sources: water is fed from an open-surface source such as canal, river, or a 
reservoir, often through a sump and intake. The water level may change over a large range and 
sediment trapping structure may be necessary. 

b) Water supply from treatment plants: treated water from treatment plants is usually supplied to a 
distribution network or a storage tower-reservoir through a pumping station. The water is clear and 
free from sediments, hence no screen and sediment settling structure is needed. 

c) Storm water pumping: storm water is full of suspended sediments and a coarse screen should be 
installed before the pump’s intake. The pumping station is used intermittently. 

d) Sewage (untreated) pumping: a sewage pump should be able to pass all solid matter through its 
system. Stagnant areas or corners must be avoided. 

e) Abstraction from borehole: the pumps used are normally less bulky (around 100-400 mm diameter), 
which could be fitted into well diameters of 150-600mm. 


17.2 Turbines 
A turbine is a device that converts the energy in water into rotating mechanical energy. 


1) Classification of Turbines 
a) Impulse Turbines: use water velocity to move the runner, rather than pressure (e.g., Pelton Wheel). 


Download free ebooks at bookboon.com 


161 


Please click the advert 


Concise Hydraulics Hydraulic Machinery 


b) Reaction Turbines: mainly use pressure rather than velocity (e.g., Francis turbine and Propeller 
turbine). 


A Pelton Wheel is a disc with buckets attached to the outside edge. The jet strikes the buckets one at a 
time, causing the wheel to spin. For a Francis turbine, water is introduced just above the runner and all 
around it and then falls through, causing it to spin. A Propeller Turbine resembles a boat propeller 
running in a tube and operates on a similar principle. Type ‘Pelton Wheel’, ‘Francis turbine’ and 

‘ Kaplan turbine’ in Wikipedia for further information. 


Selection of turbines is based on the following factors: performance (i.e. efficiency at various 
discharges), energy production, machine cost, availability, powerhouse construction cost, etc. 
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2) Hydraulics of Turbines 


The geometry of turbines is specifically shaped so that the fluid exerts a torque on the rotor in the 
direction of its rotation. The shaft power generated is available to drive electric generators. The 
hydraulics of a Pelton Wheel is used here to illustrate the application of the momentum principle in 
turbine analysis. The wheel has a series of split buckets located around its periphery. When the jet 
strikes the dividing ridge of the bucket, it is split into two parts that discharge at both sides of the 
bucket. The speed of the wheel is kept constant under varying load through use of a governor that 
actuates a mechanism that changes the setting of the nozzle. 


bk 
% 


The momentum force is 
F = pQl(V -u) cos f-(V —u)] 


The force from water to the blades is 
F,, = pAV(V -u)(1- cos f) 
P=F u= pAV(V —-u)(1-cos #)u 


There is no power developed when u=0 or when u=V. For a given turbine and jet, the maximum 
power occurs at an intermediate u that can be found by differentiating the equation above and equating 


to zero. 
< = pAV(V —2u)(1—cos f) =0 
u 


from which u=V/2. Thus the greatest hydraulic efficiency (neglecting fluid friction) occurs when the 
peripheral speed of the wheel is half of the jet velocity. 
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Power 


0.0 0.5V 1.0V 


Blade speed u 
The maximum power 


Py. = 0.25pAV°(1—cos f) 


Torque 
Tyra = FR = ReAV(V —u)(1-cos £) 


Pelton wheels are best suited (i.e., most efficient) for low flow rate and high head operations. 


The efficiency of turbines is defined by 
_ power delivered to the shaft Tao 


power taken from the water - yQh 


where T is the torque delivered to the shaft by the turbine, wis the rotative speed in radians per second, 
Q is the discharge, and h is the net head on the turbine. The efficiencies of various types of turbines 
change with load. The impulse turbine maintains high efficiency over a wide range of loads, with 
significant decrease in efficiency occurring when the load drops below about 30% of normal load. The 
efficiency of propeller turbines is very sensitive to load, and a significant drop in efficiency is 
experienced as the load falls below the normal load. In contrast, the Kaplan turbine, with its adjustable 
blades, maintains high efficiency over a wide range 
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17.3 Pump and Pipeline 


1) Performance curves 


Discharge ~ Head (Q~hp) 
The following curves can be found using a 
pump test (or obtainable from the Pi — 
manufacturers) ee eee 


Centrifugal en 5 


Discharge and head hp~Q, 7g 100) 
Discharge and power Ps ~ Q < 
Discharge and efficiency n~Q 


The discharge and efficiency curve is not 0 20 40 60 80 100 = 120—Ss«*£40 
. . . j fe) 

independent and is derived from the other two Puscuatee Weal 

curves. 


life- 
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Discharge and Head 
The discharge and head curve is the most important curve among the pump characteristics. It 
demonstrates that a pump can deliver a Discharge ~ Power (Q ~ P) 

wide range of discharges but there will be 
the changes in pressure as the discharge 
changes. The speed of at which the pump 
runs can also change the discharge and 
head. Each curve is produced under a ahaa 
specified speed and several curves should 
be obtained from the manufacturers under 
different speeds so that a whole picture 


can be obtained. In the curve, the normal 0 20 40 60 80 100 120 140 
operating condition is represented as 100% Discharge (%) 
and all other values are relative to this point. It can be seen that all pumps have higher heads at low 


discharge and lower heads when discharges go up. 


Discharge and Power 


All pumps need power to rotate their Disehacgoxutiiciency 


impellers. The amount of power needed 


depends on the speed of a pump and the i 
head and discharge required. For centrifugal 
pumps the power requirement is low when 


4 Centrifugal ie 


starting up but it rises steadily as the 


Efficiency (%) 


discharge increase. For axial pumps, the 
trend is different. There is a very large 


power demand when starting up and it will 0 


decrease as discharge increases. Mixed flow : 7 ie °° 7 a me si 
. . ] fe) 
pumps lie in between other two types and Diseiaerae) 


their power consumption is more uniformly distributed. 


Discharge and Efficiency 

It is quite important to consider the pump efficiency under its operational conditions, because of losses 
occurring in transferring fuel energy to water energy via the power unit and pump. The losses are 
caused by friction and water turbulence. The maximum range of efficiency is usually between 30 and 
80% so there is only limited range of discharges and heads over which pumps operate at maximum 


efficiency. 
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2) Initial pump selection 


Many features of a pump are summarised into one factor, i.e., specific speed, which provides a 
common base for comparing pumps. It is the speed at which a pump will deliver 1m’/s at 1.0m head 
when operating at peak efficiency and is calculated as 


sone 


Ss 3/4 
h, 


where n is rotational speed of the pump (rpm), Q is the pump discharge (m*/s), H is the pump head (m). 


Specific speed is independent of the pump size so it describes the shape of the pump rather than how 
large it is. 


Specific speed for different pumps (Kay 1998) 


Pump type | Specific speed ns Comments 

Centrifugal | 10~70 High head, low discharge 

Mixed flow | 70~170 Medium head, medium 
discharge 

Axial flow | >170 Low head, large discharge 


A user can calculate the specific speed required and check the table to find the suitable pump types. 
3) Similarity and pump model 


From dimensional analysis, the scaling law can be derived, so it is possible to experimentally 
determine the performance characteristics of one pump in the laboratory and then use these data to 
predict the corresponding characteristics for other pumps with the geometrically similar shape under 
different operating condition. The two geometrically similar pumps will have the same specific speed. 


For two pumps with the geometric ratio A = L/L, 


= AP ae , (mis the angular speed of the pump rotor) 


For the same pump, 
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2 3 


Example: A centrifugal pump having an impeller diameter of 1m is to be constructed so that it will 
supply a head rise of 200m at a discharge of 4.1m’/s of water when operating at a speed of 1200rpm. 
To study the characteristics of this pump, a 1/5 scale geometrically similar model operated at the same 
speed is to be tested in the laboratory. Determine the required model discharge and head rise based on 
the pump similarity law. 

Solution: 

The scale 4,;=5 


From 
Qugt, Qu nu = 41 _0.0328m°/s 
Oy Ay nA; 5 
From 
h ' h > 200 
=A a , Ayy = 5 “| = 7 = 8m 
hay Dey Ay \on 5 
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4) Pump selection and pipeline 


Duty Point 

A duty point refers to the point in terms of head and discharge at which a pump normally operates. 
There are many pumps on the market and there may be several that can meet the duty point or near it. 
The next is to examine the efficiency of each pump at the duty point and select the one that will operate 
near the maximum efficiency. 


The pipeline could have a few diameter options and idasimnraietfieisndy-semae 


several pipeline curves could be generated by the hy 


hydraulic loss equation. Select the diameter which 
Duty points 


could make the duty point within the maximum 


efficiency range. Several pumps and various 
an Pump curve 
Pipeline curves 


; ‘ : Sais ; (various diameters) 
analysis considering the initial capital cost and q 


running cost of the whole system should be used as Matching a pump to a pipeline 


diameters should be compared and economic 


the final decision criterion. 


Example: Water is to be pumped from a river through a 150mm diameter pipeline 950m long to an 
open storage tank with a water level 45m above the river. A pump is available and has the discharge ~ 
head performance characteristics shown below. Calculate the duty point for the pump when the friction 


factor A=0.04 (ignore local losses). 


Pump head (m) 30 50 65 80 


Discharge 2000 | 1750 | 1410 | 800 
(litre/m) 


Solution: 
Calculate the pipeline curve with A=0.04, L=950m and D=0.15m 


2 

5 ae 

D2g 
For different discharge, we can derive 
Q hf hf+45m 
2000 46.0 91.0 
1750 35.2 80.2 
1410 22.9 67.9 
800 7.14 52.4 500 1000 1500 2000 2500 


Q (l/min) 
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The duty point can be found at Q= 1400 I/s, H= 65m 


Series and parallel arrangement 

When the required flow or head is greater than a single centrifugal pump, two or more pumps could be 
combined to meet the demand. Also if the demand varies greatly with time, it is desirable to have a 
few smaller pumps that can be run in parallel instead of a single one. 


To increase the head, two or more identical pumps are operated in series. The same flow passes 
through Pump 1 and Pump 2 so that the discharge is the same as individual pumps but the head is 
doubled (in case of two pumps). The discharge~head curve for two pumps can be obtained by taking 
the curve for one pump and doubling the head for each value of discharge. Like in series electricity 
circus, the whole system breaks down when one pump stops. 


Head double for two pumps 
H 


Suction pelivery 
; pipe 
pipe : 
— 
One pump 


Two pumps in series 


To increase the discharge, two or pumps could be operated in parallel. Again the pumps must be 
identical, have individual suction pipes but share a common delivery pipe. The head is the same as for 
a single pump but the discharge is doubled. The dischare~head curve is obtained by taking the curve 
for one pump and doubling the discharge for each value of head. Usually one pump is used to provide 
low flow demand and other pumps are used to provide extra flow when the demand is increased. 


Suction 


pipes 


discharge doubles 


Delivery —> 
pipe 


One pump ie: 


Two pumps in parallel 


Download free ebooks at bookboon.com 


170 


Please click the advert 


Concise Hydraulics Hydraulic Machinery 


Cavitation and safe installation height 

When a pump is operating, it draws water from the 
sump and the water is pushed by the atmospheric 
pressure. At the sea level, the atmospheric pressure is 
about 10m of water head and would be considerably 
less in high mountains. Considering the hydraulic loss 
along the pipeline, the practical limit is below 7 meters. 
In addition, low pressure will increase the possibility of 


cavitation within the pump, which occurs when the 
liquid pressure at a given location is reduced to vapour 
pressure of the liquid. When this occurs, vapour bubbles 
form (the liquid starts to ‘boil’). Cavitation can cause a 
loss in efficiency as well as structural damage to the pump. To decide the maximum elevation for a 
pump installation, an important parameter is introduced: NPSHp (Net Positive Suction Head), which is 


a value provided by pump manufacturers as a reference for limiting suction lift under each discharge to 


ensure that a pump operates satisfactorily. 
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The actual NPSH for an installation can be calculated as NPSH,=P, - Z- V> - hr 


where P, is the atmospheric pressure (m), Z is the elevation difference between the pump and intake 
water level (m), V, is the vapour pressure of water (approximately 0.5m, depending temperature) and 
h; is the head loss in the suction pipe (approximately 0.75m, local loss is more dominant). To prevent 
cavitation, it is necessary that 


NPSH , > NPSH, 


Example: 

A centrifugal pump is to be placed above a large open water tank with discharge 14 litres/s. The 
required net positive suction head (NPSHR) is 4.6 m, as specified by the pump manufacturer. If the 
water temperature is 27°C and atmospheric pressure is 101.353 KPa, determine the maximum height Z, 
that the pump can be located above the water surface without cavitation. Ignore the friction loss, and 
the local loss is mainly caused by a special inlet filter with K,=20. The diameter of the inlet pipe is 
10cm. 


The water properties are 


Temperature (°C) 0 4 10 20 30 40 

Density (kg/m”) 999.9 | 1000.0 | 999.7 | 998. | 995. | 992.2 
2 7 

Vapour pressure 610.5 | 872.2) | 1228 2338 | 4243 | 7376 

(N/m?) 


Solution 
2 2 2 
V _ 99 L0:014/(0.05 1)] 


2g 2g 
From the table, water at 27°C, density = 996 kg/m’, vapour pressure = 3671 N/m’, 


The energy loss h, = K, =3.24m 


From NPSH,=P, - Z- V, - hy= NPSHR 
So 
101353/(996g)- Z- 3671/(996g)- 3.24 > 4.6, 


7 = 2 lem 
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Surges in pumping stations 
Water hammer can cause damages to the 


pipelines and this could also be a problem for 1. Pump stops suddenly and water 
continues to flow to create water column 


pump mains, which is the slower mass saa iia 
movement of water often resulting from the 
faster-moving water hammer shock waves. 
When a pump stops, the flow in the pipe does 


not stop immediately but continues to move 


along the pipe. Since the driving force is gone 2. Flow reverses and a sudden 


and the water will gradually slow down due to eset 
friction. As the flow moves away from the 
pump and no water enters the pipe, the empty space forms near the pump. This is called water column 
separation and pressure in the empty space drops rapidly to the vapour pressure. Then the water begins 
to flow back towards the pump gathering speed as it goes, and comes to a sudden stop as it hits the 


pump (similar to a sudden closure of a valve on a pipeline). This could burst the pipe and the pump. 


Possible solutions are: 

1) Stop pumps slowly: This will avoid the water column separation and minimise the surge. Diesel 
pumps are more suitable for this solution since they stop gradually after the fuel is cut off, while 
electric pumps usually have sudden stop when the electricity is switched off. 


2) Use a non-return valve: it will allow water to pass in one direction only and protect the pump (not 
the pipeline) when water flows back. 


3) Use an air vessel: This is similar to a surge 
tank. When a pump stops and pressure starts 
to drop, water flows from a pressured tank into 


the pipeline to fill the void and stop the water 


Water flows in and out of vessel 
column from separating, and water is allowed 


to cushion water hammer 
to flow back into the tank when the flow is 
reversed. The water will oscillate back and 
forth until it eventually stops through friction (much like a car shock absorbers). This device is 

expensive and should be used only for large scale pump stations where serious water hammers are 


expected. 
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Questions 17 


Hydraulic Machinery 


1. Each day 600m? of water is pumped by 10m to a tank above the ground. Calculate the amount of 
power required to do this (assume the system is 100% efficient and there is no energy loss in the 
pipeline). 

(Answer: about 0.7 kW) 


2. A Hydraulic power plant takes in 30m*/s of water through its turbine and discharges it at V=2m/s at 
atmospheric pressure. The head loss in the turbine and conduit system is 20 meters. Estimate the 
power extracted by the turbine. 

(Answer: 23.5 MW) 


= Elevation 150 m 


Wales Turbine 


-— 2m/s — Elev. 50m 
—— 


3. A power plant is to be built on a stream capable of utilising a flow up to 1 m’/s with a gross head of 
40m at 78 percent overall efficiency. Assume the stream has enough flow to maintain the plant 
working at full capacity continuously and the power is valued at £0.05 per kilowatt-hour. What annual 
revenue could be expected from this plant? 

(Answer: £134,000) 


4. A Pelton Wheel is to be selected to drive a generator at 600rpm. The water jet is 75mm in diameter 
and has a velocity of 100m/s. With the blade angle at 170°, the ratio of vane speed to initial jet speed 
at 0.47, and neglecting losses, determine 1) diameter of wheel to centerline of buckets (vanes); 2) 
power developed. (Note P = pAV(V —u)(1—cos f)u ) 


(Answer: 1.497m, 2.18 MW) 


Download free ebooks at bookboon.com 


174 


Please click the advert 


Concise Hydraulics Hydraulic Machinery 


Solutions 17 
Hydraulic Machinery 


1. Convert discharge to m’/s 
Q =600/ 24/3600 =0.007m’ /s 
P= pgHQ =1000x9.81x10x 0.007 = 686W = 0.7kW 


2. Set up the energy equation between the reservoir surface and the turbine outlet 


V2 V2 
yzt+ta=4t+t+h+H 
28 2g 
2 


2 
so 100=— +h, +H 
2g 


With h=20m, so H =100- 20 -0.2 = 79.8 m 


Power extracted by turbine will be 
P= pgQH = 1000 x 9.81 x 30 x 79.8 = 23.5MW 


6007 'A'g SweIsAS WIy] JER] @ 


iB \ 
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3. Power output P = npgQH/1000 kw 
= 0.78x1000x9.8 1140/1000 
=306.1 kW 
Annual revenue=Px365x24x0.05=£134.1x10° 
So the annual revenue from this plant will be 134 thousand pounds. 


4. 1) The peripheral speed of the wheelis =u=047x100m/s=47m/s 


u _ Xe) | _ A7 
a 60=600rpm D 5007 x 60 107 1497m 
2) The poweris P= pAV(V —u)(1—cos f£)u 


2 
sg. P= 1000 =~°= 100(100-47)(1- cos170°)47 


= 2183154W = 2.18MW 
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Appendix: Further Reading Resources 


The following resources are highly recommended if you want to further explore the interested topics in 
hydraulics. This is not an exhaustive list and will be updated regularly in the future. You are welcome 
to recommend useful books/web sites that are not on the list. 


Biswas, A.K., 1970, History of Hydrology, North-Holland Publishing Company 


Chadwick, A., Morfett, J. and Borthwick, 2004, Hydraulics in Civil and Environmental Engineering, 
4" Edition, Spon Press 


Chanson, H., 1999, The Hydraulics of Open Channel Flow, Butterworth-Heinemann, Oxford 

Chow, V.T.,1959, Open Channel Hydraulics, McGraw-Hill College 

Evett, J.B. and Liu, C., 1989, 2500 solved problems in Fluid Mechanics & Hydraulics, McGraw Hill 
French, R.H., 1985, Open-Channel Hydraulics, McGraw-Hill 

IAHR, 2008, Media Library, http://www.iahrmedialibrary.net/ 

Kay, M., 1998, Practical Hydraulics, E&FN Spon 

Liggett, J.A. and Caughey, D.A., 1999, Fluid Mechanics, ASCE Press 


Munson, B.R., Yong, D.F. and Okiishi, T.H., 2002, Fundamentals of Fluid Mechanics, John Wiley and 
Sons. 


Novak, P., Moffat, A, Nalluri, C. and Narayanan, 1996, Hydraulic Structures, Second Edition, E & FN 
Spon 


Wikipedia, 2008, ‘Hydraulics’, ‘Hydrostatics’, ... http://en.wikipedia.org/wiki/Hydraulics 
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